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ABSTRACT 


Previous investigations of longshore currents have included various simplifving as- 
sumptions and restrictions including a planar beach, a steadv and depth uniform flow. 
spatially-invariant bed shear stress and turbulent momentum exchange. and the exclu- 
sion of surface wind stress. These assumptions are quantitatively investigated bv calcu- 
lating the relative importance of each term in the longshore momentum balance with an 
emphasis on the relative importance of wind forcing across the barred nearshore. Wind 
and wave forcing of longshore currents across a barred beach are examined using both 
a numerical model and field measurements. A local momentum balance was measured 
at various locations across the surf zone during the SUPERDUCK experiment held at 
fee rym. Corps of Engineers, CERC Field Research Facility. Duck. N.C.. U.S.A. 
in October 1986. A moveable sled was instrumented with pressure. current, and wind 
sensors to measure the various terms in the longshore momentum equation. Stabilitv- 
dependent atmospheric drag coefficients for the surf zone are determined from wind 
stress Measurements acquired just beyond the surf zone and wind speed measurements 
acquired from an anemomieter atop the 9 m sled mast. Breaking waves Were visually 
identifted and electronically marked on the data tapes. Mean current velocities at three 
levels were measured across the barred surf zone to examine the effects of a bar on the 
memtical structure of the tlow. A modthed longshore current model is used to study the 
effects of barred topography. The predicted longshore current profiles are compared with 
field data and bed shear stress coefficients are determined from model-fitting. Significant 
spatial variability of longshore current [O(1.0 ms~!)] is observed across the bar with 
depth uniform longshore flow. Wund force is found to be a first-order term along with 
wave force under certain wind and wave conditions. Spatially-dependent bed shear 
stress coefficients relative to the bar are determined from a local momentum balance. 
Mean bed shear stress coefficients were 0.004 + 0.0010, 0.002 + 0.0006, and 0.001 + 
0.0003 for offshore the bar, on top and immediately before the bar, and in the trough. 
Mean bed shear stress coefficients determined from model-fitting and from the momen- 
tum balance are in excellent agreement. Longshore current velocities are modeled within 
meee 7o Of Observed velocities across the barred beach at SUPERDUCK with a mean 


least square error of 7-15%. 
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I INTRODUCTION 


Previous investigations of longshore currents have utilized many simplifving re- 
strictions and assumptions. Major assumptions and restrictions include a planar beach, 
a steady and depth uniform longshore current, spatially-invariant bed shear stress and 
turbulent mixing, and the exclusion of surface wind stress. These assumptions are often 
too restrictive to permit realistic results. This research quantitatively investigates these 
assumptions by measuring the relative importance of each term in the longshore mo- 
mentum balance equation (see equation (1)). 

Modern theories for mean longshore currents are based on a longshore momentum 
balance (Bowen, 1969; Thornton, 1970; and Longuet-Higgins, 1970a, 1970b) with a ra- 
diation stress approach (Longuet-Higgins and Stewart, 1964), For a plane-parallel 


beach, the complete longshore momentum equation 1s: 
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Where a right-handed Cartesian systein 1s adopted with x increasing offshore (see 

Figure 1), and 

term | = temporal change of the time-averaged depth-integrated mean mo- 
mentum per unit area in the longshore direction due to both steady (MM) 
and unsteady flow ) 

term 2 = cross-shore gradient of the wave-induced longshore momentum flux 
(also called the wave-induced radiation stress gradient), 

term 3 = bottom shear stress in the longshore direction, 

term 4 = surface wind stress in the longshore direction, and 

term 5 = cross-shore gradient of the longshore momentum flux due to turbu- 


lence. 


Most investigators (Longuet-Higgins, 1970a; Jonsson et al., 1974; Kraus and Sasaki, 
1979; Thornton and Guza, 1986; Losada et al., 1987; Baum and Basco, 1987; among 


others) apply the previously mentioned assumptions and restrictions and then reduce the 
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Figure 1. Section and plan definition for the SUPERDUCK experiment (mod)- 
fied from Birkemeier and Dalrvmple, 1975). 
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equation for longshore momentum to a simphfhied balance between the wave-induced 
radiation stress gradient, the bottom shear stress, and the turbulence-induced radiation 
stress gradient. The most restrictive of these assumptions 1s that of a planar beach, since 
barred beaches are more common than planar beaches (Wright et al., 1986). Temporal 
variability of longshore currents, often exceeding + 100% of the mean flow, have been 
measured in the field (Dette and Fuhrboter, 1974; Wood and Meadows, 1975; Meadows, 
1976). The depiction of spatially-invariant bed shear stress and turbulent mixing across 
the dynamic surf zone, especially over a barred beach, may be unrealistic. Additionally, 
recent investigations addressing combined currents and waves indicate that the bed shear 
stress coefficient, which is directly proportional to bed shear stress itself, may be 1.5 to 
3.5 times larger in magnitude (Grant et al., 1984) than the invariant values universally 
used by longshore current investigators. 

Longshore current forcing is studied utilizing the data sets acquired during the 
SUPERDUCK experiment. This experiment was held at the U.S. Army Corps of Engi- 
neers, Waterwavs Experiment Station’s Coastal Engineering Research Center (CERC), 
Field Research Facility (FRF) in Duck, NC in September and October 1986. A unique 
instrumentation system designed specifically for this experiment was utilized to measure 
a local momentum balance at various locations across the surf zone of this single-barred 
beach under various wave, current, and wind conditions. Instruments were mounted on 
a mobile sled to measure the various terms in the momentum balance of equation (1). 
Two triangular, differential pressure sensor slope arrays Were spaced approximately 2 m 
apart. Coupled with an absolute pressure reference, the slope arrays can be used to 
measure the gradient of the wave-induced momentum flux. Marsh-McBirnev 
electromagnetic (EM) current meters located at the centroid of the slope arravs can be 
used to determine both mean momentum as well as the total, or combined, momentum 
flux due to the wave-induced and turbulent components. The difference between the 
current meter determination of the total momentum flux (S,,) and the slope array de- 
termination of the wave-induced momentum flux (Sy) yields a measure of the depth- 
integrated Reynolds stress, or turbulent component of the momentum flux. The surface 
Stress can be inferred from the wind velocity measured using a wind anemometer 
mounted atop the 9 m sled mast. 

Three methods for calculating the bed shear stress coefficient are attempted -- 
measuring the momentum balance, model-fitting to field data, and the mean profile 


method. By operating in areas where a plane-parallel beach assumption is valid and by 


precisely orienting dynamuc measurements perpendicular to the local bottoin contours, 
all the terms contributing to the longshore momentum can be measured. Having meas- 
ured all the terms in the balance, the unspecified bed shear stress coefficient is then de- 
termined. Theoretical mean longshore current distributions generated by a numerical 
model are compared with field-measured longshore current speeds. The bed shear stress 
coefficients in the numerical model are then adjusted until the modeled and observed 
current distributions coincide in a least square sense. Preliminary indications from an 
on-going surf Zone investigation suggest that the mean longshore Velocity profile in the 
vertical appears to be logarithmic (personal communication, R. T. Guza, 1987). Hence, 
the bed shear stress coefficient might also be determined by the profile method of 
boundary layer physics. However, logarithmic profiles were not found in the SUPER- 
DUCK data. 

The sled was designed to make transects anywhere along the shoreline and as far 
offshore as the mast height would allow. Sled mobility was provided by U.S. Army 
Corps of Engineers equipment. The sled was tethered by a 300 m chain to a heavy-duty 
forkhft on the beach. The sled was towed offshore by the Coastal Research Amphibious 
Buggy (CRAB) to the desired location beyond the surf zone and then detached from the 
Sled. The forkhft dragged the sled onshore along its pre-planned transect of data col- 
lection locations. 

A previous difficulty with directly measuring momentum flux in the field was the 
precise orientation required of the instrumentation (Thornton and Guza, 195))) Tiga 
SUPERDUCK experiment, precise orientation of the sled (and attached slope arrays) 
was achieved to within 0.5° using laser ranging to triangulate on two prisms niounted 
2.4m apart on a mast spreader above the sled. A daily calibration of current meter 
orientation was accomphshed by iniually positioning the sled beyond the surf zone for 
the first data run of a surf zone transect. Based on conservation of radiation stress 
outside the surf zone, the current meters Were numerically oriented so that their respec- 
tive radiation stress values were identical. This orientation was then held constant for 
the remaining sled locations along that transect. This calibration procedure is essential 
to obtaining realistic physical results from the data and is a distinguishing feature of the 
research. 

A primary improvement between the SUPERDUCK experiment and) previa 


large-scale nearshore experiments was the extensive bathymetric data. Daily 
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bathymetric profiles of the entire experiment area were accurately measured with the 
CRAB and the laser ranging system. 

Pertinent literature is reviewed in the following chapter. A numerical model 1s de- 
veloped in Chapter II] to compare wind and wave forces for various wind speeds and 
directions. The SUPERDUCK field experiment ts described in Chapter 1V. Data anal- 
wsis of the field experiment is extensive and therefore 1s addressed in separate chapters. 
Preliminary data analvsis is described in Chapter V and includes an overview of the ob- 
served mean current structure over the bar. Since surface wind stess is of major impor- 
tance to this research and few wind stress studies have been conducted tn the surf zone, 
Chapter VI is entirely devoted to the parameterization of the surf zone surface wind 
stress for this field experiment. Calculation of the longshore momentum balance, de- 
termination of the relative contributions of each term in the balance, and calculation of 
spatially-variable bed shear stress coefficients are addressed in Chapter VII. A longshore 
current model for planar beaches 1s modified to account for barred topography and 1s 
applied to the data. The model's longshore current results are compared to observations 
and bed shear stress coefficients are determined by model-fitting im Chapter VIII. 
Summary and conclusions are presented in Chapter IX. Because of the large quantity 
of data presented graphically for each dav of data analysis, onlv a single typical dav 
(October 16) 1s illustrated within the text. Data analysis figures for other days are ref- 
erenced in the text by figure number but are collectively grouped in Appendix A. 

ive data obtained from this field experiment represents the first quantitative data 
set. including directional wave. surface wind stress. bathymetry. and breaking wave 
measurements, for the study of longshore current over barred beaches. In most previous 
studies. wind forcing of longshore currents has not been considered due to the absence 
of quantitative data. This research quantitatively determines the role of wind forcing in 
driving longshore currents from both field measurements and a simple numerical model. 
Wind forcing, as well as a variable bed shear stress coefficient, can represent measurable 
and reasonably significant terms in the momentum balance equation, and therefore must 


be taken into account in any quantification of longshore currents. 


Ca 


HI. LITERATURE REVIEW 


A. INTRODUCTION 

One of the earliest works on coastal hydrodynamics and processes was that of 
Johnson (1919). He qualitatively discussed wave set-up, longshore currents, and rip 
currents (undertow) as hydraulic currents due to waves. He utilized incident wave angle 
as his criterion for rip currents (waves approaching perpendicular to shore) and for 
longshore currents (waves approaching oblique to shore). 

In World War I], amphibious landing craft had difficulty landing their troops in the 
presence of a longshore current. The current often represented a greater hazard than the 
breakers themselves by causing the craft to broach inside the breaker zone. [iiemiames 
Navy's concern for this effect prompted Putnam et al. (1949) to conduct one of the first 
applied research studies on longshore currents. 

Since that time, various theories have been postulated to account for the basic 
physics of longshore currents and their predicted velocity distribution across the surf 
zone. These theories include conservation of energy, conservation of mass, empirical 


formulation, and conservation of momentum . 


B. RADIATION STRESS 

Modern longshore current theorv is based on conservation of momentum. One of 
Putnam et al.’s (1949) approaches was the apphcation of conservation of momentum. 
The momentum flux of an obliquely incident wave was averaged over a wave period to 
give a mean momentum flux entering the surf zone. The longshore component of the 
change (due to wave breaking) in this momentuin flux provided a driving force for 
longshore currents. They assumed a steady-state longshore current, and therefore the 
driving force was balanced by a bottom frictional force. 

The momentum flux approach was formahzed and the physics clarified when the 
concept of radiation stress was introduced. Longuet-Higgins and Stewart (1964) defined 
radiation stress as the “excess flow of momentum due to the presence of the waves.” 
Elaborations on the concept of radiation stress are contained in Longuet-Higgins and 
Stewart (1960, 1961, 1962, 1963, 1964), Bowen (1969), Longuet-Higgins (1970a, 1970b, 
1972a, 19120), and Vhomtond(t 3 10" 

Theory was further formahzed with the derivation of time-averaged and depth- 


integrated conservation of mass and conservation of momentum which separates the 


mean and fluctuating momentum contributions (Phillips, 1977; Mei, 1983). Thus ap- 
proach has been used extensively in the study of wave set-up and set-down, horizontal 
longshore current profiles, rip currents, and surf beat generation. This approach forms 
the foundation of modern longshore current theory, as Well as the foundation for the 
momentum balance of equation (1) on which this dissertation is based. 

When there is an oblique wave approach to the shoreline, the longshore component 
of radiation stress (or longshore momentum flux due to unsteady motion) is given by 


(Longuet-Higgins, 1970a): 


n 
Syx = [. pu'v'dz (2) 
Miere 7 = the sea surface. 

h = water depth as measured from the still water line (SWL). 

p = water density, 

uv = fluctuating (x. y) components of the flow velocity, and 


the overbar denotes a time-averaging operator. 


Mie Oasic Concept is that the longshore current velocity 1s related to the change in in- 
coming wave momentum, which in turn relates to the radiation stress. This change of 
momentum flux in the cross-shore direction is proportional to the rate of wave energy 
dissipation, where dissipation is due primarily to surface breaking but can include bot- 
tom, and internal (turbulent) dissipation. By using linear wave theory, radiation stress 
for monochromatic waves can be related to the incident wave angle (a) at the measure- 


ment location by: 


eae? Sina cos a (3) 


where E is total wave energy density, a is the wave incident angle, or angle between the 
wave crest and a line parallel to the shoreline, and n is the ratio of wave group velocity 


to wave phase velocity expressed in terms of wave number (k) as: 


a 2h 
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Radiation stress has also been apphed to random waves (Collins, 1972; Battjes, 1972a; 
Thornton and Guza, 1986) and non-linear wave theorv (James, 1974a, 1974b) to solve 
for the longshore current. 

The conservation laws of mass, horizontal momentum flux, and radiation stress by 
Phillips (1977) and Mei (1983) are both widely apphed and similar. Both develop the 
equations for mean water motion by integrating over depth and averaging over a wave 
period. They differ in their approach to the mean flow due to unsteady motion (1. 


waves). Separating the mean and fluctuating horizontal velocity components (u;) as: 


u =U; + i i=1,2, (5) 
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Thus he incorporates the mean flow due to unsteady motion in the mean velocity term, 


U,, resulting in a mean momentum per unit area defined as: 


Awe — M,; = pUuy + h) (S) 

Phillips (1977) utilizes equation (5) and (6), however he does not include the mean 
flow due to unsteady motion in the term U;. He evaluates the time-averaged depth- 
integrated current flow by partitioning the integration of u (equation 7) into two sepa- 


rate regions: 
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Where the first term on the right-hand side is zero by definition. This results in a mean 


momentum) per unit area of 


CO 


n 
M,=pU,(7t+h) + p| uwdz i= 1,2 (10) 
0 
Phillips approach is adopted for this research. 
Since the momentum balance of equation (1) represents the basis for this research, 
each remaining term -- surface wind stress, bottom shear stress, turbulent momentum 


flux, and temporal variability -- will now be examined. 


C. SURFACE WIND STRESS 

Surface wind stress is a shear stress imparted to the water column when the wind 
flows across the water surface. For simplicity, surface wind stress has been historicallv 
neglected as a forcing function for longshore currents. Although no physically-based 
model has specifically evaluated the wind stress contribution to longshore currents, there 
is empirical evidence suggesting wind stress can be important, and therefore should be 
included as a driving force for longshore currents. Nummedal and Finley (1978) used a 
stepwise multiple regression procedure utilizing surf zone environmental parameters to 
explain the variability of longshore current velocities. Although this procedure was 
empirical and not theoretical, the longshore current component of the wind velocity was 
found to account for more of the observed current variance than any single parameter, 
or combination of parameters, descriptive of the breakers. 

Fox and Davis (1973) were able to relate longshore currents with barometric pres- 
sure at Lake Michigan. Although wind is certainly related to barometric pressure, thev 
did not expand their work to directly deal with wind. Komar (1976) mentioned that 
winds blowing in the longshore direction may contribute to the observed flow. 

Birkemeier and Dalrymple (1975) are one of the few nearshore modelers who did 
incorporate a wind forcing term (albeit, an underweighted term) in their longshore cur- 
rent model. Notwithstanding the different spatial scales involved, it 1s noted that lake 
(e.¢., Sheng et al., 1978) and offshore modelers do consider wind stress whereas 
nearshore modelers do not. Nearshore modelers either assume that wave stress domi- 
nates over surface wind stress or that no wind is present. 

Hubertz (1986) calculated the spectra of the mean longshore wind component and 
mean longshore current for a four-month period at the FRF in Duck, N.C. Longshore 
Current spectral peaks were present at periods of 4.6 days, 23.8 hours, and 12.5 hours. 


Longshore wind component spectral peaks were present at 5.9 days and 24.4 hours. He 


related the two longest period current peaks to wind forcing by synoptic weather systems 
and by sea breeze conditions, respectively. The 12.5 hours peak was attnbuted tomtiaes 
There are several methods for calculating surface wind stress on the water sumieas 


These methods are: 
e the direct method (also called the covariance or eddy correlation method), 
e the profile method, 
e the dissipation method, and 


e the drag coefficient method (also called the bulk aerodynamic method). 


Each method has its own inherent advantages, disadvantages, and accuracy. The drag 
coefficient method was selected for this research based on equipment availability, the 
rough environment and platform involved (surf zone measurements by a mobile sled), 
and the availability of wind stress measurements acquired by other field investigators at 
the SUPERDUGK cxpemmen: 

The drag coefficient method formulation for surface wind stress (Tavlor, 1916) 1n the 


longshore direction 1s: 
Y , , 
= Pica (11) 


where 7t% = time-averaged surface wind stress vector in the longshore (¥) direction 
(hereafter referred to as surfacemina stress). 
Po = atmospheric density, 
C, = an atmospheric drag coefficient (generally of O( 107) for a height of 10 
m above the water surface) for the specific height above the water surface 


at Which W is measured, 


W = the total wind speed, and 
it’, = the longshore component of the wind velocity. 


The absolute value is used in equation (11) to retain directionality. 

A large quantity of research has been directed at determining values for C, for land 
and the open ocean. Garratt (1977) provides an excellent review for investigations of 
C, for tand and the open ocean prior to 1977. More recent investigations include Wu 
(1980, 1982), Safare (1984), Geernaert and Katsaros (1986), Geernaert et alee 
1987a, 1987b), and Hsu (1986a, 1986b, 1987). Unfortunately, few investigations have 


dealt with determining C, near the land-sea interface and in the surf zone. C, magnitudes 
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peo 7 10 to 2410) were detennined for nearly calin conditions outside tlie 
surf zone of a Chukchi Sea barrier island (Walters, 1975) and Cj; magnitudes of O( 
2.0 x 10-3 to 7.0 x 10-3) for the upper swash zone of a Florida beach (lsu. 1970). 
Walter’s (1975) investigation was conducted under adiabatic and onshore wind condi- 
tions. Whereas Hsu's (1970) investigation included adiabatic as well as non-adiabatic 
conditions with onshore wind conditions. 

fst ( 1973) determined suri zone wind stress values for onshore winds of 0.5 to 8.5 
ms—! under adiabatic conditions for a two-week period at Santa Rosa Island, FL. Cy, 
magnitudes of O(1.0 x 10-2 to 1.5 x 107) are obtained upon conversion of his surface 
wind stress values. No breaker wave height values were indicated; however, 1t 1s assumed 
that they were small due to the experiment’s location (northwest Florida coast), month 
(June), wind speeds, and a photograph aceompanying the published results. 

Vugts and Cannemeyer (1981) caleulated C, for non-adiabatie and onshore wind 
conditions from field measurements at a Danish beaeh. They instrumented two 20 m 
Masts with multi-level cup anemometers and temperature sensors and placed one at the 
high-tide line and one directly inland at a distance of 127 m. Thev concluded from ap- 
proximately 100 half-hour profiles diat the drag coefficient during unstable conditions 
over the surf zone is the same as over the ocean. Thev caleulated C, magnitudes of O( 
ex 10-* to 1.5 x 10-4) for wind speeds at 10 m height from 3-10 ms7!. 

Geernacrt ect al. (1987) measured Wind stress using a sonic anemometer from a re- 
searcl: platform in the North Sea. They found an empirical relation between the neutral 
drag coefficient (C,,,) and the ratio c,/u., where c, is the dominant long wave phase speed 
and uw, 1s the surface wind friction velocity. They apphed this relationship to other data 
sets taken in varving water depths. They found increased Cy, values for shallower depths 
(Figure 2). As wind speed increased to 25 ms7!, the C,,, value for 10 mm depth was in- 


Seedsed bY a fuctor of 2 (to O(3.0 x 10-4)) over that for deep water. 


D. BOTTOM SHEAR STRESS 

Putnam and Johson (1949) conducted one of the earliest investigations of nearshore 
bottom stress or bottom friction. They assumed sinusoidal oscillatory bottom motion, 
an invariant friction coefficient, negligible perpendicular flow at the bottom due to 
percolation, a planar bottom slope, and a constant water density. Based on their inves- 


tigations, bottom shear stress (7°) may be represented by: 


ee |i, (12) 


1] 


10 15 20 25 


WIND SPEED (m/sec) 





Ficure 2. Calculated dependence of Cz, on the 10 m height wind speed for five 


depths (from Geerneart et al., 1987b), 


where Gy the dimensionless bed shear stress coefficient due to waves and currents 
and generally of O( 107+), and 
u, = resultant vector sum of the mean longshore current velocity (Vv) andamge 


instantaneous wave orbital speed near the bottom (see Figure 3). 


As with wind stress, the absolute value is used in equation (12) to retain directionality. 
This equation has essentially formed the basis for modern bottom friction studies. In 
general, later investigators have worked toward its modification as well as solving for 
empirical values for cy. 

One can decompose u and v into a steady or mean velocity (U, V) and an unsteady 
or fluctuating velocity (u’, v’). The unsteady velocity may be further separated into 


wave-induced (u, ¥) and turbulence-induced (u”, v”) components (Mei, 1983): 
pe ae ee (13) 


ve Vty = Vey + yp" (14) 


Incorporating the above decomposition, the bottom shear stress is represented as a 
vector oscillating in both direction and magnitude due to the interaction of a quasi- 
steady current with superimposed wave motion. From conservation of mass flux and the 
assumption of straight-and-parallel bathymetric contours, the mean cross-shore current 


velocity (U) is assumed Zero. uw, 1s then expressed as: 
- + ‘lie alaee 
ett +(V+tuy =u ty Sut i} 


mere ul. = uy, COSG, 


ety, Sila, 
u,, = Wave orbital speed just above the bottom and immediately outside the 


bottom boundary laver, and 


A 
i, &FR unit vectors in the x and v directions, respectively. 


Bottom shear stress (t?) can be written as: 
7b Bee Be 
mot, I + TJ (16) 


and from equations (12) and (15): 


E pc, u, lu + pcy| uly (17) 


The longshore component of time-averaged bottom shear stress can then be written as: 
b 2, ab 
ty = perluly = pfu + )2y¥ (18) 


The time-averaging procedure must be performed after evaluating the velocity compo- 
nents. 

The bed shear stress coefficient is a function of bottom roughness and flow intensity 
near the bed. Bottom roughness is a function of both sand characteristics and bed form, 
where density and grain size define sand characteristics and bed form is described by the 
existence of ripples, which are in turn, a function of the flow field. Hence, the bed shear 
Stress coefficient, bottom roughness, and flow intensity are all interrelated. 

Often in geophysical research, laboratory results complement field measurements 
and assist in building a particular data base as well as increasing our knowledge of the 
phenomena being studied. Unfortunately, very few of the large number of laboratory 


studies of cy are appropriate for extrapolation to the field. The flow field in nature will 
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almost always be in the rough turbulent flow regime (Jonsson, 1966). Since the 
Reynolds number is the appropriate dynamical similitude scale parameter, large 
Reynolds numbers are required to be in this flow regrme. Hence, only laboratory studies 
performed in large wave flumes or oscillating water tunnels, where longer wave periods 
can be obtained, appear appropriate. 

Numerous investigators have solved for cy under waves (Table 1) with results rang- 
ing from O(107! to 10~) for different beaches. Most have sought an invariant value for 
cy based on a given set of conditions (e.g., bottom roughness, beach slope, or for only 
one particular beach), without regard for any spatial or temporal variability. Vitale 
(1979) states that a common engineering practice 1s to assume cy = 0.01. Johns and 
Jefferson (1980) showed the variation of bottom stress within a wave cycle in a numer- 
ical analvsis. They found increasing cy with increasing bottom roughness and decreasing 
wave period. 

There are four methods of bottom shear stress, and subsequent cy , determination: 
-- momentum balance, model-fitting, the profile method, and direct shear measurements. 
The first three are indirect methods and are addressed below. The direct shear meas- 
urement method involves directly measuring the shear force, such as with pressure plates 
fixed on the bottom or measuring the Revnolds stress at the bed. This has not been 
done 1n a surf zone and 1s therefore onutted from further discussion. 

1. Momentum Balance 

Forcing terms and assumptions are specified (often based on data availability) 
to form an appropriate longshore momentum balance equation. The bottom stress term 
and its accompanying bed shear stress coeffictent are then determined as residual values. 
The strength of this method depends on how complete the momentum balance ts and 
how realistic are the restrictions and assumptions. Thornton and Guza (1981) utihzed 
this method at Torrey Pines Beach, CA (a near-planar, unbarred beach) in connection 
with the Nearshore Sediment Transport Study (NSTS) experiments. Assuming steadv- 
State conditions, neglecting wind forcing, and operating in an area of stratght-and- 
parallel contours, they assumed a simple balance where the cross-shore gradient of the 
total longshore momentum flux (both wave-induced and turbulent components) was 
balanced by the mean longshore bottom shear stress. For simplicity, they neglected the 
temporal variation term. They found considerable temporal and spatial variation of cy 
with a mean value of 0.01 + 0.01 . However they stated that the results were to be used 


with caution due to the difficulties in directly measuring the radiation stress (S,.,) and its 
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Table 1. 


BY VARIOUS INVESTIGATORS 


Investigator 


Bretschneider and Reid (1954) 
Kish! (1954) 

Iwagaki and Tsuchiya! (1966) 
Twagaki and Kakinuma (1967) 
Jonsson (1966) 

Hasselman and Collins (1968) 
Longuet-Higgins (1970a) 
James? (1972) 

Jonsson et al. (1974) 
Kamphius (1975) 

Komar- (1975, 1976) 

Hiumties (1976) 

Liu and Dalrymple (1978) 
Hsiao and Shemdin (1978) 
Madsen et al.? (1978) 

Kraus and Sasaki (1979) 
Vitale (1979) 

Johns and Jefferson (1980) 
Grant et al. (1984) 

Wind and Vreugdenhil (1986) 
Guza et al. (1986b) 

Thornton and Guza (1986) 


of 

0.01 

0.03 - 0.04 
0.01 - 0.04 
OC Pee 
0.01 

0.015 

0.01 

O-0025 70:00) 
0.05 10510) 
0.00420. 250 


OO Sas 
0.0026 + 0.0006 


0.04 
0.002 - 0.100 
0.0] 
O.O11 - 0.024 


0.003 - 0.250 
0.0096 - 0.00319 
0.005 - 0.011 
0.0015 

0.008 

0.006 + 0.0007 


BED SHEAR STRESS COEFFICIENTS PROPOSED AND UTILIZED 


Comments 


Various data sets 


f(bottom roughness) 
non-linear eqn. 


various data sets 


«if required, the investigator's friction factor was converted to a bed shear Stress 
coeilicient freni¢-={imction factor 4). 


1 from Sonu (1975) 
2 from Basco (1982) 


gradients. These difficulties were due to angularity errors caused by lack of instrument 
resolution, installation misalignment, and errors associated with choosing the correct 
alongshore direction. All of these errors prevented accurate directional resolution of 
waves with small angles of incidence (< 3°). Resolution of these difficulties was ad- 
dressed in the SUPERDUCK experiment by utilizing a precise surveying SvStemameume 
selecting davs for analvsis when the mean wave incident angle was large (> 8°). In the 
following determination of cy by the momentum balance method, both the temporal 


variability and wind stress terms in equation (1) are included. 


2. Model-fitting 

Theoretical horizontal mean longshore current distributions generated by a nu- 
merical model are compared with field-measured longshore current speeds. The cy in the 
numerical model is then adjusted so that the two profiles agree. The resulting numerical 
model values for c, can then be assumed to apply for the field measurements. Selecting 
the appropriate physics for the numerical model 1s obviously important for this method. 
Thornton and Guza (1986) used this method with a narrow-band, random-wave trans- 
formation model for data collected in another NSTS experiment at Leadbetter Beach, 
CA. A hierarchy of cy solutions were derived for linear and non-linear formulations for 
bottom shear stress, and with and without turbulent momentum exchange. Thornton 
and Guza (1986) found a least-squares “best fit” value of cp= 0.006 + 0.0007 using the 
full non-linear. and presumably more correct, formulation for the bed shear stress. 

3. Profile Method 

Miiiesprolemethodwe.c., sec Roll, 1965; Schlichting, 1968) has been utilized in 
boundary laver physics for vears to determine friction or shear velocities and roughness 
heights. It assumes a constant stress layer with a mean logarithmic Vertical flow struc- 
ture and requires at least two simultaneous speed values at different heights above the 
same location. 

Grant et al. (1984) conducted high-quality near-bottom boundary layer meas- 
urements at a continental shelf location in 90 m water depth. Even at this substantial 
depth, they found the near-bottom flow field to be composed of oscillatory currents due 
to swell superimposed on mean currents, a situation similar to Figure 3. The mean, 
near-bottom (<2 m) velocity profiles were found to be logarithmic. Using the profile 
method, mean stress values were 3 to 7 times larger than expected from predictions of 
mean current alone using a typical smooth-bottom drag coefficient, and 1.5 to 3.5 tunes 
larger than expected for predictions using a drag coefficient based on the observed rough 
bottom. This significant increase in mean stress magnitudes was attributed to the pres- 
ence of the swell wave oscillatory currents. These results may apply to the surf zone, if 
the surf zone’s entire vertical extent of chaotic flow may be likened to a fully developed. 
neutral, turbulent boundary laver in the longshore direction. 

Analogous to surface wind stress (e.g., see Roll, 1965; Schlichting, 1968), the 


longshore bottom shear stress may be written as: 


t= pr (19) 
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where 1+ 1s the longshore shear or friction velocity associated with the mean flow. 


Combining equations (18) and (19): 


> 2 
Vs Vx 


(a (20) 
(uty) Gp 2 ey 


The time-averaged velocity profile in a neutral turbulent boundary layer over a 


topographically simple bottom 1s logarithmic and may be expressed as: 


V(2)=(—)In | (21) 


where xk is von Karman’s constant and and z, 1s the roughness height. By simultaneously 
measuring current speed at two different heights (V(z,) and V(z,)), two equations in the 
form of (21) may be written. The two unknowns, ¥ and z,, are first solved, and then c, 
is obtained from equation (20). Grant et al. (1984) employed a similar profile method 
in their deep water (90 m depth) investigation where friction velocity was measured in a 
downstream direction and u and v were measured in the streamwise and cross-stream 


directions respecuvelx. 


E. TURBULENT MOMENTUM FLUX 

The horizontal transfer of momentum due to turbulent miximg processes is not well 
understood. When the classical analytical solutions for the horizontal profile of the 
longshore current Were initially solved using monochromatic breaking waves (Bowen, 
1969; Longuet-Higgins, 1970a, 1970b; Thornton, 1970), a velocity maximum and dis- 
continuity was observed at the breaker line, with zero longshore velocity at the shoreline 
as well as bevond the breaker line. This occurred since there were no driving forces 
outside the sharply defined breaker line to generate a longshore current and there was 
maximum forcing occurring at the breaker line. However, field and laboratory data 
showed a smooth transition of the longshore current from outside the surf zone to a 
maximum inside the breaker line. This observation was attributed to turbulent mixing 
(lateral shear stress). Analytical solutions for longshore currents on a planar beach were 
proposed with the assumptions of isotropic turbulence, a planar beach, monochromatic 
waves, and the parameterization of a lateral shear stress in terms of an eddy viscosity 
coefficient multiplied by a horizontal velocity gradient. Values for the eddy viscosity 


coefficient were obtained by model-fitting. 


ss 


Whereas the bed shear stress coefficient is a physical quantity related to bottom 
roughness and flow intensity, the eddy viscosity coefficient is a function of the flow field 
and is a parameterization used to remove the longshore current discontinuity at the 
breaker line. An alternative approach to a monochromatic wave model is to describe 
breaking wave heights in terms of a probability distribution. This results in a more re- 
alistic region of wave breaking rather than a distinct breaker line. This spreading out 
of the breaker line results in a smooth horizontal profile of the longshore current. thus 
eliminating the need to include an eddy viscosity (Collins, 1970; Battjes, 1972a). 

Thornton and Guza (1986) utilized a random wave model and model-fitting to de- 
termine values of the bed shear stress coefficient and eddv viscosity coefficient for both 
a linear and non-linear expression for the bed shear stress. Eddy viscosity coefficients ( 
Bereere represented by an adjustable coefficient, N, in the range of O< N < 0.016, 
expressed by (Longuet-Higgins, 1970a, 1970b): 


ae 
2 


mee | xi(gh)2 x> 0, (22) 


where g is the acceleration due to gravity. Optimal cyand N coefficients were determined 
by least-squares fitting between calculated and measured longshore current speeds 
shoreward of the mean breaker line. Using a non-linear bed shear stress formulation. a 
cp Value of 0.006 was determined. Eddy viscosity was found not to be important for the 
near planar topography of their field experiment. 

The last term in equation (1) represents an expression for horizontal diffusion or 
lateral mixing of turbulent momentum in terms of a depth-integrated Reynolds stress. 
The longshore component of momentum flux due to unsteady motion (equation |) 1s the 
depth-integrated covariance between the fluctuating velocity components of u and Vv, 
measured by the current meters. By partitioning the unsteady motion as that due to 
waves and turbulence (equations 13 and 14), the excess momentum flux (equation 2) can 


be separated into a wave-induced (S,) and turbulent component (S",,): 


= 7 on q a ted 
ee yx t 5 yx = ik puvdz + i, pu'v dz (23) 


It has been assumed for this separation that the characteristic time scales of the two 
fluctuations are vastly different (ie., u; and u”, are uncorrelated and 4, u", = OQ). inte- 
grated viscous stresses are negligible, and that the horizontal scales of depth and mean 


velocity are comparable. The turbulent component of the excess momentum flux due 
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to unsteady motion may be obtained by comparing the total radiation stress (Sie) 
measured by a current meter to the wave-induced radiation stress (S,,) measured by a 


slope array. 


F. TEMPORAL VARIABILITY 

Most previous investigators have assumed steady-state conditions for longshore 
momentum. However, significant temporal variations in longshore current velocity 
measurements have been observed (Putnam et al., 1949; Inman and Quinn. 1951; Inman 
and Bagnold, 1963; Dette and Fuhrboter, 1974; Dette, 1974; Huntley and Bowen, 1975; 
Wood and Meadows. 1975; Meadows, 1976, 1978: Guza and Ihornton, 197S-)Gigaaes 
al., 1986b, and others). The temporal variability of mean currents (o\4,/67) will be in- 
cluded in the analysis where Af, is the trme-averaged and depth-integrated mass flux per 
unit area, or momentum per unit area, 1n the longshore direction due to steady (M,) and 


unsteady (/,’) flow: 


) 


M,= 3h + My = pV7th)+ | _ pv'dz (24) 


Inman and Bagnold (1963) ascribed temporal variability of longshore currents to the 
wave groupness. Meadows (1976) concluded that the total observed longshore current 


velocity was composed of three components: 
fag (3 Ps a a 
v=] steady are wave 7 y long (23) 


where Voeady 18 the steady longshore component, V,,,,. 1s the coincident wave period 
fluctuating component, and I,,, 1s the long period (as related to surf beat or low mode 
edge waves) fluctuating component. He defined the steadv state component as the ve- 
locity component whose period of variation was greater than 200 seconds. The selection 
of 200 seconds was not explained. Guza and Thornton (1978) concluded that a temporal 
averaging period for mean longshore currents was not well-understood. Various re- 
searchers utilize different averaging trmes to determine the steady state component, 


based more on record length or data analysis, than on rationale. 


G. VERTICAL DISTRIBUTION OF MEAN LONGSHORE FLOW 
The equations by Phillips (1977) and Mer (1983) mntegrate the water motion over 
depth. For simplicity, the longshore current is assumed uniform from the surface to the 


bottom. No longshore current theories account for a depth-dependent longshore cur- 


rent. Theoretical work on the depth-dependency of longshore currents (Stive and de 
Vriend, 1987) and cross-shore currents (Svendsen et al., 1987) has only recently been 
undertaken. Most laboratory and field investigations have not measured the depth- 
dependence of longshore currents, and as a result, there is little field data on vertical 
current structure. 

Greenwood and Sherman (1983) measured vertical structure of the longshore cur- 
rent at the Wendake Beach, Ontario field experiment. They calculated mean flows for 
three elevations as well as a depth-averaged mean longshore current velocity. Their 
maximum mean variability was onlv 10% of the depth-averaged flow. Thus their data, 
as well as Meadows’ (1978) and Hubertz’s (1986) data, suggests that modeling longshore 
currents as uniform over depth 1s appropriate. Although they assumed a logarithmic 
boundary laver velocity profile in order to calculate a boundary roughness length (z,), 
their logarithmic plots of only two points are inconclusive regarding the existence of a 
logarithmic vertical velocity profile in the surf zone. However, their linear plots of ver- 
tical profiles show a “consistent deformation of the velocity field strongly suggestive of 
the influence of bed friction,” and hence a logarithmic profile within at least 0.5 - 1.0 m 
of the bed remains a possibility for their data. 

Wright et al. (1986) obtained field data of vertical profiles from 0.1 to 1.0 m above 
the bottom of longshore currents in a bar trough. Although the longshore velocities 
Were on the same order of magnitude as the wave orbital velocities, the bottom bound- 
ary laver was dominated by the waves acting 1n the shore-normal direction, orthogonal 
to the longshore currents. The velocity profiles were not logarithmic and it appeared 
that they were strongly influenced by wave-generated vortices from the bed ripples 
present. 

Theoretical investigations by Lundgren (1972), Smith (1977), and Grant and 
Madsen (1979), laboratory studies by Bakker and Van Dorn (1978) and Grant and 
Madsen (1982), and field measurements by Grant et al. (1984) assume that the mean 
velocity profile is comprised of an inner region (close to the bed which is directly affected 
by wave-induced turbulence) and an outer region above it, and that the outer region has 
a logarithmic profile which is also affected by wave-induced turbulence. Preliminary 
indications from: an on-going surf zone field investigation suggest that the mean 
longshore velocity profile in the vertical might be logarithmic (personal communication, 
R.T. Guza, 1987). The existence of a logarithmic profile would allow bed shear stress 


coefficients to be determined by the profile method. 


H. BARRED BEACHES 

Although barred beaches are more common (Wright et al., 1986), far less is known 
about surf zone currents on barred beaches than on planar beaches. Allender and 
Ditmars (1981) found little spatial variability in the horizontal longshore current profile 
across the entire surf zone of a steep, barred beach during storm conditions on Lake 
Michigan. No significant cross-shore flow or vertical variability of current flow was 
observed. Greenwood and Sherman (1983, 1986) studied surf zone currents in Lake 
Huron off a steep, multi-barred beach at Georgian Bay. They found significant cross- 
shore variability in the longshore current, and up to a four-fold variation during storm 
conditions. These two contrasting conclusions are indicative of the embryonic state of 
knowledge concerning the surf zone dynamics of barred beaches. The SUPERDUCK 
experiment was conducted on a barred beach. Data analvsis of the surf zone currents 


will be used to improve the understanding of surf zone dynainics at a barred beach. 
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II. WIND FORCING OF LONGSHORE CURRENTS FROM A 
FINITE-DEPTH SPECTRAL WAVE MODEL APPROACH 


To assess the importance of the wind contribution to longshore currents, wave and 
wind forcing terms in the longshore momentum equation are integrated across the surf 
zone in a simple model and compared. Wave and wind forcing are coupled by assuming 
the waves driving the longshore currents are generated by the local wind. The TMA 


wind-wave shallow water wave spectrum is used as input to the radiation stress. 


A. THE TMA SPECTRUM 

The TMA spectrum (Bouws et al., 1985) was developed from data obtained from 
two extensive shallow-water wave measurement programs. MARSEN (Marine Remote 
Sensing Experiment at the North Sea) and ARSLOE (Atlantic Ocean Remote Sensing 
Land-Ocean Experiment), and data from a prolonged storm near lightship TEXEL in the 
southern North Sea. The MARSEN and ARSLOE sites were both on the continental! 
shelf with depths less than 40 m. The combined data set was named TMA, after the 3 
sub-data sets, and includes water depths from 6 - 42 m, wind speeds up to 30 ms“, 
bottom composition from fine to coarse sands, and beach slopes ranging from !:150 to 
nearly flat. The uniqueness of the TMA data set to this research is that the ARSLOE 
data, which comprises two-thirds of the nearly 3000 wind-wave spectra, were measured 
emeae CERC Field Research Facility in Duck, N.C. 

Bouws et al. (1985) proposed, and then verified, that a finite-depth wind-wave 
spectrum could be approached through application of the similarity principles bv 
Kitaigorodski et al. (1975). The TMA spectrum is an extension of the self-similar 
JONSWAP (Joint North Sea Wave Project) spectral shape (Hasselmann et al., 1976) to 
include the effects of finite water depth embodied in the Kitaigorodskii et al. (1975) 


finite-depth transformation factor, D,(w, h): 
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and 
0. = Phillips constant, 
Ww = 2rf, 
f = frequency, 
Le = peak spectral frequency, 
vy = peak enhancement factor. and 
Oo = ¢, = peak spectral width on low frequency side of f,,. 


= o, = peak spectral width on high frequency side of /,. 
Eps Se/m.&) 1s the Pierson-Moskowitz (1964) spectral shape for fully-developed deepwater 
waves, and O,(ff, yy. Gg,0,) is the JONSWAP spectral shape function. 

Thompson and Vincent (1983) developed an approximation for @,(w,/) using linear 


theory and redefining the two independent variables (@ . h) as a single variable (w,): 


> wo? Wp, = | 
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However, the approximation (equation 30) is incorrect for w, > 2 and therefore is mod- 


ified here as: 
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@(w,) approaches a value of one in deep water and a value of zero as the depth de- 


creases (Figure 4). 





Figure 4. Kitaigorodskil et al. (1975) finite-depth spectral transformation function, 


®,.(w.h) (from Bouws et al.. 1985), 


Bouws et al. (1985) developed the following empirical parameter expressions for in- 


put into the TMA data set: 


90,0078 K’*” (33) 
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Where /,) = wind speed atz = 10 mand 
ke = wave number for waves at peak frequency. 


These expressions are valid for all water depths in the TMA data set. Thus, the TMA 


spectrum requires specification of #¥,, , h, and /,. 


B. NUMERICAL APPROACH 

To compare wind and wave forcing for various conditions using finite-depth spectra, 
it 1s assumed that: 

1} Wave conditions are solely the result of a local wind (i.e. there is no swell 
present), and therefore the wind direction (@) and mean incident wave direction 
are identical, and 

2) Sea conditions are fullv arisen (1.e., the sea is not duration-limited). 

Comparisons are made for inputs of gentle and steep planar beach slopes (tan f = .01 
and .10), wind speeds (10-30 ms~!), all wind directions, and peak spectral frequencies 
conimensurate with wind speed. 

Surface wind stress is determined from equation (11). C, 1s determined from meas- 
urements taken by Geernaert et al. (1987b) in open waters of the North Sea for a depth 
of 10 m (Figure 2). These C, values are increased bv 33% to account for surl Zoi 
roughness. The 33% increase will be addressed later in the Chapter VI. Peak frequency 
is determined for the Pierson-Moskowitz spectrum (Carter, 1982): 


ley 


In = (37) 





Peak frequency 1s assumed not to shift during shoaling (Hughes, 1984). 

After input of tan f and H¥’,5, peak wavenumber (K,,) 1s calculated using the linear 
theory relationship starting at an initial water depth of 20 m. The TMA spectrum 1s 
determined from equation (26). The RMS wave height is obtained by integrating the 


TMA spectrum: 


0.5 7 
(HimstMa = 8) Ema ‘ (38) 
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where only sea waves (maximium frequency of 0.5 Hz) are considered. 
Thornton and Guza (1982, 1983) found an envelope curve relating //,.,, to depth in 


the inner surf zone: 


fe. = 0.42 h (39) 


The f7,,,, determined by equation (38), which includes the effects of variable depth, is 
checked against equation (39) to ascertain if wave breaking due to shoaling had com- 
menced. If the ratio H,,,,/h is less than required for surf zone wave breaking, water 
depth is reduced and a new H,,,. determined. This process for calculating breaker depth 


iterates until equation (39) is satisfied. Surf zone width (x,) 1s then determined from: 


h 
tan £ 





xp = (40) 

The radiation stress gradient (@S,, / Ox) and surface wind stress (ty) are both inte- 
grated with respect to x, to allow comparison of a wave force (S,,) and a wind force 
(ti xp), with units of force per unit longshore length. The integration is confined to x, 
since the radiation stress gradient 1s conserved outside the surf zone where there is no 
wave breaking. Momentum exchange between the wind-driven longshore current out- 
side x, and the wind- and wave-driven longshore current inside x, 1s assumed negligible. 


Using equation (2), wave force is expressed as: 
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where the energy is expressed in terms of variance, which is: 


0.5 

oT MAxn(f = yO n(f) Af (42) 
0.0 

and n(f) is expressed in a frequency-dependent form. The factor of 1/2 is an adjustment 

factor for broad-banded wind-driven waves which incorporates the effect of short- 

crestedness on the ratio of radiation stress to the average energy content of the waves 

(Battjes, 1972b), 


Wind force 1s constant over the surf zone and is determined from equation (11): 


Xx 
tx, = [3 dx = p, Cyl WIV sin 6 x, (43) 


C. RESULTS AND DISCUSSION 

Resultant spectra, wind force, and wave force for both a gentle beach slope and a 
steep beach slope are illustrated as cases A through F in Figure 5 and Figure 6. H,,,., 
h, far Wind speed, Cy, and x, for each of the six cases are listed in Table 2. In¢@reaame 
wind speed effect for a gentle beach slope (tan f = 0.01) is depicted in cases A, B, and 
C, and for a steep beach slope (tan f = 0.10) in cases D, E, and F. Wind and wave 
direction are relative to normal beach incidence. The three spectra previously described 
(Pierson-Moskowitz, JONSWAP, and TMA) are included to illustrate the effect of 
shoaling on the wind-wave spectra, as Well as the expected shift to lower peak spectral 
frequency as wind speed increases. The //,,,. and depth (h) at commencement of 


breaking also increase with increasing wind speed. 


Table 2. CASE PARAMETERS FOR WIND VERSUS WAVE FORCE COM- 
PARISONS — (see Figure 5 and Figure 6). 
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The wind force is directly proportional to the width of the surf zone (x,), which 1s a 
function of tan f. As a beach steepens, surf zone width narrows and thus there is a re- 
duced area for surface wind stress to act. Wave force is independent of tan f. 


Wind direction dependency is inherent in equations (41) and (43) where: 


S.J cc sin@ cos@ = —~ sin(26) (44) 
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Figure 5. 
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Comparison of wind and wave force for a gentle beach slope (tan fp = 


0.01) with associated wave spectra. 
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Figure 6. 





.00 





0.25 
Freq (Hz) 





0.50 


Force (sim ”) Force (Ism”*) 


Force (Jsm’*) 





é —180 -—30 0 30 180 
(Gh = 
10d 
10 
oa . 
an 

: —180 -—30 0 90 180 
10 


ar 


Nay Kan ee ( | Ds Cede Ta 
~ bade 


©), Oh 
ul il dina 


©) 


— = 
= 
=e. 
- = 
ewe @e 


l 
Co 
oO 

l 
wo 
oO 
wo 
o) 
Oo 
>) 


Comparison of wind and wave force for a steep beach slope (tan B = 0.10) 


With associated wave spectra. 


30 


ty oc sin 6 (45) 


Thus the maximum force occurs at + 90° (alongshore) for wind and + 45° for waves. 
Wind force increases more rapidly than wave force for increasing 6. When winds are 
blowing offshore (i.e. absolute wind directions > 90° ), waves do not form because of 
fetch limitations and thus only wind force exists, creating a weak, but sorally wind-driven 
longshore current. 

The dependency of wind and wave force on wind speed Is not as readilv apparent 
because of the interrelationship of wind speed with the various TMA parameterizations 
(equations 33 to 37). The absolute force of both wind and waves increases as wind speed 
increases. The increase in wave force 1s due to the increase 1n spectral energy. The in- 
crease in wind force 1s through the wind speed and the increase in surf zone width caused 
by larger waves (generated from the TMA spectrum) breaking farther offshore. Thus 
the wind force’s dependency on wind speed ts indirectly dependent on the spectral wave 
energy. 

Wave and wind force are compared for increasing wind speed and constant wave 
incident angle (Figure 7). Since wave force 1s independent of # . onlv wind force 1s 11- 
lustrated for both gentle (solid line) and steep (dashed) beach slopes. The amount of 
increase for both wind and wave force over the wind speed range of 5 ms~! to 30 ms7! 
is nearly four orders of magnitude. Wand force 1s seen to increase more rapidly than 
wave force for increasing wind speed and constant incident wave angle. Varving the in- 
cident wave angle results in the same conclusion; therefore only @ equal to 20° 1s illus- 
trated. 

The more rapid rate of increase for wind force 1s quantified in terms of a wind force 
to wave force ratio: 
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The ratio of forces is depicted in Table 3 (tan 8 = 0.01) and Table 4 (tan f = 0.10) for 
the quadrant of wind direction from 0 — 90°. The ratio is symmetric for wind directions 
of 0 to —90° and is a function of slope, wind speed, and wind direction. For steep beach 
slopes, wind force is much less than wave force and thus may be neglected. However, 
for gentle slopes, wind force increases as the wind becomes more atongshore in direction. 


For a given wind direction, an increase of 10 ms7! in wind speed results in a 50% in- 
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Figure 7. Comparison of wind and wave force for mncreasing wind speed at constant 


wind direction (0 = 20°). 
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crease in the ratio. Although this ratio for gentle beach slope is generally less than unity 
for all wind speed and direction combinations, wind force is never less than an order of 
magnitude compared to wave force. Wind to wave force ratio versus beach slope for a 
typical mid-latitude cold frontal passage is depicted in Figure 8. Wind and wave frontal 
conditions are defined as a wind speed of 15 ms! and @ = 20°. The wind to wave force 
ratio increases with decreasing beach slope. The ratio is approximatelv 10% for a beach 
slope of 0.04 and approaches 40-50% for beach slopes less than 0.01. Therefore wind 
force should be considered in studying longshore currents during strong winds on gentle 
beaches and during offshore winds. In Chapter VII, wind to wave force ratios from the 


field experiment will be compared to the results in Table 3 and Table 4. 


Table 3. RATIO OF WIND FORCE TO WAVE FORCE FOR A GENTLE 
BEACH SLOPE (tan B = 0.01) 


Wind Direction (relative to normal incidence) 


e JNF - ratio is infinite since there is no wave force 


® waves are assumed to be locally generated and to approach the shoreline from 
the wind direction 
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Table 4. RATIO OF WIND FORCE TO WAVE FORCE FOR A STEEP BEACH 
SLOPE (tan ff = 0.10) 


Wind Direction (relative to normal incidence) 


e INF - ratio is infinite since there 1s no wave force 


@ waves are assumed to be locally generated and to approach the shoreline from 
the wind direction 
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Figure 8. Wind to wave force ratio versus beach slope for a typical mid-latitude cold 


frontal passage. The wind and wave incident angle (@) 1s relative to 


beach normal. 
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IV. FLELD EXPERIVEE Na 


A. INTRODUCTION 

To develop an improved understanding of coastal processes (currents, waves, 
sediment transport, and nearshore geomorphology) under a wide variety of conditions, 
a one month-long field experiment was held in the autumn of 1986 at the U.S. Army 
Corps of Engineers, Waterways Experiment Station’s Coastal Engineering Research 
Center (CERC), Field Research Facility (FRF) at Duck, N.C. Duck 15 on thei 
Ocean side of the North Carolina Outer Banks and is located approximately 10 miles 
north of Kitty Hawk, N.C. This major field experiment was entitled SUPERDUCK and 


included nearly 20 scientific investigations. 


B. SITE DESCRIPTION 

The FRF site (Figure 9) is located along a 100-km unbroken stretch of shoreline 
of an extensive barrier island formation known as North Carolina’s “Outer Banks.” The 
barrier island upon which the FRF site is located is bordered by the Atlantic Ocean to 
the east and Currituck Sound to the west. Except for a few fishing piers, therevareume 
major coastal structures or littoral barriers along the entire reach of shoreline. The site 
has a tidal range of 0.5 to 2.0 mand regular offshore bathymetry free of features which 
may affect wave climate. The beach is characterized by a mean foreshore slope of 1:10, 
a single longshore linear bar svstem which can become three-dimensional, and a mean 
slope of 1:100 offshore of the bar. Shoreline excursion due to tidal fluctuations is typi- 
cally 10 m. Mean sediment grain size is approximately 0.5 mm. 

FRF personnel support includes a staff of nine (scientists, technicians, divers, and 
heavy equipment and computer operators). Extensive equipment support includes a 
LARC-V amphibious vehicle, a Coastal Research Amphibious Buggy (CRAB), auto- 
mated data processing (ADP) equipment, a heavy-duty forklift, four-wheel drive and 
all-terrain Vehicles (ATV's), laboratories, electronic workshops, a maintenance garage, 
and high-resolution Zeiss Elta-2 electronic surveying systems. 

Several planning workshops were held months prior to the experiment to synergize 
the efforts of all the principal investigators. The FRF staff and SUPERDUCK partic- 
ipants jointly installed more than 70 electronic sensors in the nearshore environment for 
the experiment. The overall set-up of the various instruments at the FRF site is illus- 


trated in Figure 10. The longshore current meter array was approximately 500 m in 
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length and 25 m offshore. Its primary purpose was for the investigation of infragravity 
waves, however its data will also be used in this research to investigate temporal and 
spatial variability of longshore momentum flux and to provide an extra data station for 
sled transects of the surf zone. Unfortunately, the cross-shore array of current meters 
and pressure sensors was not functional for the experiment. The third cross-shore sen- 
sors Were sonar bottom profilers. A high-resolution linear pressure sensor array in 8 m 
depth provided offshore directional wave data. Iligh resolution bathymetry data were 
collected on a nearly daily basis between profile lines 165 and 275 by laser positioning 
the CRAB with the Zeiss surveving system as it methodically traversed the area along 
the profile lines. Example products of the surveying effort are included in Figure 11 and 
Figure 12. A major difference between this and previous large-scale nearshore exper- 


iments was the extensive bathymetric data acquired using the CRAB. 


C. EQUIPMENT DESIGN AND OPERATION 

Only a brief description of the sled and its instrumentation is provided here. 
Martens and Thornton (1987) describe in detail the design, instrument calibration, and 
specification of the sled system. The all-aluminum sled (Figure 13) was constructed by 
the CERC Field Research Facility staff in Duck, N.C., based on a Naval Posteradiae 
School design. The sled foundation 1s two runners 3 m in length on which was placed 
a3 mx 3m aluminum sled frame with a 9 m aluminum sailboat mast. The frame was 
approximately one quarter of a meter above the sled runners. All! the instrumentation 
was then attached to this frame. Significant additional weight was placed on the runners 
to provide stability and prevent dragging in the robust surf zone. The sled was designed 
with a low profile to minimize flow disturbance. 

The differential pressure slope array concept was patterned after Bodge and Dean 
(1984) with some major modifications. Two-inch O.D. stainless steel threaded pipe and 
pipe fittings were assembled 1n the shape of a capital “H” to form two adjacent triangular 
slope arrays with a common vertex at the center (see Figure 13). The five differentia! 
pressure ports were located at the four ends of the “I” and one at the center. The pipe 
cavity was filled with a premium grade vacuum pump oil, similar to mineral oil, which 
had a low viscosity coefficient and was assumed incompressible. Four Sensym model 
LX06002D monolithic differential pressure transducers were installed at the ends of the 
“T”, just inboard of the diaphragms which were exposed to the sea. These differential 
pressure transducers provide an output voltage proportional to the applied differential 


pressure With a response time of 0.1 seconds. The pressure at the center diaphragm was 
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Three-dimensional bathymetry example, 16 October 1986. 


Figure 11. 


tical lines represent locations of cross-shore and longshore instruments. 


40 





w “JONULSIO 


250 


CUS Wale oe: 


150 


5 


CO 
t= 
C) 
C) 
(D 
w~-4 
4 
aie 
ate 
= © 
ac 
om 
pF a 
x & 
Cc) 
i) 
OQ 
‘aa 
(1) 
QO 
= 
Te) 


Two-dimensional bathymetry example, 16 October 1986. 


Figure 12. 
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Figure 13. Side and top view of the SUPERDUCK sled. 
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hydraulically transnutted to the positive port of each differential pressure sensor, while 
the sea pressure at that location was sensed by the negative port of the sensor. A 
Paroscientific model $020 digiquartz pressure sensor was installed adjacent to the center 
diaphragm to measure the absolute pressure. 

Current velocities were measured using three Marsh-McBirney bi-directional 
electromagnetic current meters with spherical (4-cm diameter) probes mounted on the 
sled. These current meters are ideally suited for surf zone studies by providing accurate 
measurements in free stream turbulence as well as combined steady/oscillatory currents 
(Guza et al., 1986a). Both variable vertical stacking and slope array comparisons were 
considered in positioning the current meters on the sled frame. Two current meters were 
mounted at the centroids of the two triangular-shaped slope arravs at elevations 0.7 and 
1.0 m above the bed, with the third current meter midway between at an elevation of 1.5 
m. During sled deployment, the current meters Were always positioned on the “up- 
current” side of the sled so as to avoid flow contamination bv the sled mast. 

The current meters’ directional electrodes were aligned with the longitudinal axis of 
the sled frame bv eve and are considered to be within | — 2° of the true axis. During 
data analysis, a computational rotation is done for true alignment based on conservation 
of momentum flux outside the surf zone. 

A single-point R.M. Young model 05103 wind monitor was mounted atop the mast 
to provide “in-situ” measurement of vector wind in the surf zone with simultaneous wave 
momentuin measurements. 

All equipment were battery-powered by 12 sealed lead-acid batteries housed in a 
Watertight canister mounted on the sled. The analog data were PCM (pulse code mod- 
ulated) encoded and then telemetered ashore by a system described by Lowe et al. 
(1972). The data were then decoded and recorded on 9-track magnetic tape. Data ac- 
@arsition tines for all the various SUPERDUCK experiments were referenced to a pre- 
cise Eastern Standard Time (EST). 

Visual observations of waves breaking on the sled mast were marked by an electrical 
trigger which registered a pulse on the data tapes. Although extremely tedious, this ef- 
fort allows determination of the occurrence of wave dissipation events. 

isiceise Orientation of the sled, with respect to true north and the SUPERDUCK 
experiment area, was obtained with the use of the FRF’s Zeiss Elta-2 electronic survey- 
ing system. Two reflective prisms (5-cm diameter) were mounted at each end of a 2.4 


m longitudinal mast spreader located approximately 7 m above the sled. From triangu- 
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lation. the relative positioning of these two prisms and their absolute positioning within 
the FRF coordinate system provided the required precise orientation (Within 0) sagen 
current meters and slope array. Relative three-dimensional positions between various 
points on the sled and the prisms were obtained with separate Zeiss shootings prior to 
and after the sled’s deplovment in the SUPERDUCK experiment. 

Prior to selecting an area of the surf zone to transect, the previous dav S batiiviietms 
profiles and contours were studied and changes in morphology noted. The geographic 
area With the highest degree of straight-and-parallel isobaths was selected for each day's 
operations. The runs were made during daylight and scheduled for 1-2 hours before and 
after a tidal change, so as to mininuze the tidal effect and variation in mean shoreline 
location. This timing also coincided with the recording of all other SUPERDUCK data, 
which was accomplished within + 2 hours of a tidal change on a 24-hour basis. The sled 
was initially towed offshore by the FRF’s CRAB to a position just bevond the breaker 
zone for the first run of a surf zone transect. Subsequent runs were made at tegen 
of maximum breaking. on top the nearshore bar, and in the nearshore trough. Lach run 
was at least 35-minute in duration. The sted was brought inshore for each subsequent 
run of a transect by towing the sled with a chain tether attached to a four-wheel drive 
forklift located on the beach. Zeiss Elta-2 laser surveving shots of the sled’s prisms were 
taken at the beginning and end of each run, with many runs having additional surveying 
shots during the run itself. 

In addition to the sled’s anemometer, simultaneous meteorological data of air tem- 
perature, sea surface temperature, atmospheric pressure, and wind velocity were re- 
corded by other SUPERDUCK investigators at the end of the pier and near the FRF 


laboratory. 
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V. PRELIMINARY DATA ANALYSIS 


Data analysis of the SUPERDUCK field experiment is addressed in the next four 
chapters. Preliminary data analvsis 1s addressed in this chapter and includes data 
screening for data gaps and signal interference, instrument orientation, three- 
dimensional depiction of the mean current structure over the bar, and a discussion of 
homogeneity and stationarity for the experiment. Since surface wind stress 1s a major 
concern and few wind stress studies have been conducted in the surf zone, Chapter VI 
is entirely devoted to the parameterization of surf zone surface wind stress. The 
longshore momentum balance 1s addressed in Chapter VIII. The relative importance 
of each momentum term is calculated. Spatially-variable cys are determined as residuals 
from the momentum balance. A longshore current model for planar beaches is modified 
to account for barred topography and its results are compared with observations in 
Chapter VIII. 


A. INITIAL DATA SCREENING 

The analog data were PCM decoded and then digitized at § Hz. Time series plots 
for all 16 channels of information were produced for the entire data collection period. 
These records were visually perused for data gaps and signal interference. Basic statistics 
(mean, Variance. skewness and kurtosis) were calculated and formated in tabular form. 
Temporal plots of two-minute means and variances were produced. This information 
was used to eliminate bad data collection periods and data channels. Data collection 
achieved a 95.5% successful collection rate with over 40 hours of data collected over a 
10-day period. Nearly all of the eliminated data were identified at the beginning or end 
of a data run and therefore temporal continuity was maintained. 

The wind direction portion of the R.M. Young wind monitor was not calibrated 
before its unexpected destruction by wave forcing on 18 October. Therefore the mean 
wind direction measured at the sled 1s unavailable. The mean wind direction measured 
at the end of the FRF pier (SethuRaman et al., 1987) was used for this research. Mean 
wind speed after 1§ October was also taken from the end-of-pier data. 

Occasionally, one or two current meters would be partially or totally above water 
during a sled data collection period, such as when the sled was atop the bar at low tide. 
These periods were identified from visual inspection of the time series, supplemented bv 


field logbook entries. These data were not utilized in the radiation stress calculations. 


B. SLED INTER-ORIENTATION. 

A portable Zeiss laser ranging system was used to precisely measure the sled 
instrumentation positions relative to the sled. A hand-held reflective prism was placed 
at a designated measurement point on the sled (e.g., the right rear corner of the sled 
frame top) and that position designated as a three-dimensional reference point from 
which all other sled points of interest (prisms, current meter probes, slope array 
diaphragms, etc.) were measured. A complete laser measurement of the sled was con- 
ducted on 26 September 1986 and 20 October 1986. Sled modifications, such as instru- 
ment adjustment, mast replacement, etc., were noted in the field logbook so as to enable 
selection of the proper survey dav for subsequent calculations. The Zeiss measurements 
were supplemented with engineer tape measurements. 

The sled instruments Were aligned to a three-dimensional local cartesian coordinate 
system referenced to the top of the sled frame bv a rectangular coordinate rotation: 
Xnew = Xoig COS O + Hoig Sin 5 
Vnew = —*oig $11 8 + Hoig COS O 7) 
where @ is the angle of rotation. A table of absolute three-dimensional distances be- 
tween all points of interest was then generated. These distances were repeatedly used 
throughout the analvsis. 

During field operations, three-dimensional instrument coordinates, with respect to 
the FRF coordinate system were obtained from the sled prism three-dimensional coor- 
dinates by: 


e relating the mast spreader. slope array, and current meter system longitudinal and 
vertical axes to the sled frame longitudinal and vertical axes, 


e determining instrument pitch and yaw angles every 64 seconds from linear in- 
terpolation of the two-to-five survey fixes for a sled run collection period. All data 
Were then adjusted every 64 seconds for the current sled pitch and yaw angles, 


* adjusting current meter and anemometer heights above the bottom every 64 sec- 
onds for sled pitch and vaw angles, and 


¢ correcting the distance between current meters for the sled’s pitch and yaw angles. 


Since all data were taken from a mobile sled rather than from fixed-bed instrumen- 
tation, it Was imperative that the data were accurately translated into a coordinate sys- 
tem relative to the local bathymetric contours. Use of the Zeiss laser ranging svstem to 
accurately measure the sled points of interest, to fix the sled in the experiment area, and 


to provide comprehensive bathymetric data coverage, enabled this translation. 
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C. BEACH ANGLE ORIENTATION. 

Location coordinates obtained from Zeiss-sighting of the sled’s prisms were con- 
verted into sled pitch and yaw angles so as to correct the sensor data with respect to the 
horizontal and vertical and the FRF coordinate system. Pitch and yaw angles were 
calculated every 64 seconds from linear interpolation between the two most current Zeiss 
sightings to account for any sled movement during the data run. The instrument array 
orientation was then numerically rotated relative to the measured bathymethy. The se- 
lection of 64 seconds coincided with the ensemble averaging time for the spectral calcu- 
lations. 

Sled data were acquired from 11-18 October and 20-21 October 1986. Locations of 
sled data runs within the SUPERDUCK experiment site are illustrated in Figure 14. 
Incomplete survey sightings of the sled prisms due to salt incrustation on the prism lens 
and Zeiss operator error prevented determination of precise sled location on 13, 14 and 
20 October. Therefore these three davs of data were not used. 

A complete bathymetric survey was not done on 11, 17 and 21 October due to high 
surf conditions or equipment availability. For sled collection davs when no coincident 
bathymetry was surveyed, the previous or next day’s bathymetry was appropriately se- 
lected based on minimizing any storm-induced morphological change (see Table 5). 

A mean bathmetric contour direction was determined by applying a least-squares 
regression fit to the survey data offshore and located within approximately 100 - 150 m 
of the nearshore bar. This is the area of maximum bathymetric gradient with minimal 
slope farther offshore. The bathymetry inside the bar was more complex than offshore 
the bar, however the shoreward distance over which it could affect wave action was small 
O(30-50 m). The most influential bathymetry was assumed to be that located on and 
offshore the bar. 

The longshore length scale over which this regression fit was applied was weighted 
(2.3 - 1/3), with emphasis on the portion of the beach which favored the approaching 
waves (e.g., for waves approaching from the northeast, the northern offshore contours 
were favored by a factor of 2). The total distance assessed was varied from 100 to 250 
m (or to maximum extent of precision bathymetry) with resultant mean orientation 
varying less than 0.5° between various distances considered. The mean orientation of 
the four “best-fit” straight lines (corresponding to 0.5 m contour intervals) was chosen 
as the mean bathymetric contour. The beach orientation angle was the angle between 


this mean contour and the coordinate system’s y-axis. The beach orientation angle was 
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Figure I4, Sled data run locations. Open circles are sled run locations and closed 


circles are fixed-bed current meter locations. 
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then assigned as a constant for that transect and used in conjunction with instrument 


orientation which is described later. 


Table 5. BATHYMETRY SELECTION AND BEACH ORIENTATION AN- 
GLES 


Sled data dav Bathvmetrv dav Profile! Beach orientation angle (°) 2 
250 ., 
250 
255 
197 
195 
195 
170 


! Profile locations are shown in Figure 10. 

2 Beach orientation angle is measured between the FRF coordinate svstem’s 
y-axis and the mean bathymetric contour from on top and 
offshore of the bar. The angle was measured over a longshore length 
scale Varving between 100-250 m. 





A combined angle (@,) was used to convert from the yaw angle measured by the 
Zeiss to the correct rotation angle upon which to coordinate transform the u,v velocities. 
This combined angle was the sum of the sled yaw angle, the beach orientation angle, the 
angle between the current meter longitudinal axis and the sled centerline, and the angle 
between the mast spreader longitudinal axis and the sled centerline. This combined an- 
gle (@.) was used in equation (47), restated here with velocities substituted for coordinate 
locations, to perform the coordinate transformation: 

Bee = Unig COS 0. + Voig Sin 6, (48) 
me — —i74 $1 8. + Voig COS O- 

Fixed-bed current meter velocities were also converted to the FRF coordinate sys- 
tem using equation (48) where @, 1s field-measured orientation angles with respect to the 


FRF coordinate system. 


D. CURRENT METER ORIENTATION PRIOR TO A TRANSECT. 
Since the radiation stress gradient is conserved outside the surf zone, any difference 
in measured S,,(/) between the current meters on the sled is considered to be due to er- 


rors in current meter alignment with respect to each other and to the sled, current meter 
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calibration, or to imprecise orthogonal orientation of the dvnanuc measurements with 
respect to the assumed straight-and-parallel bathymetric contours. The contours outside 
the bar were nearly straight-and-parallel (see Figure 12 and Appendix B). Therefore, 
current meter alignment or gain was considered tn error. Orientation of the current 
meters’ axes with the sled’s axes was considered to bemwiiiine eee 

To correct for the radiation stress difference between two current meters outside the 
surf zone, one current meter was assumed to be true in both gain and alignment. The 
other current meter was then rotated numerically using equation (48), or adjusted in 
gain, until the radiation stress gradient was minimized. 

Two current meters and two adjustment criteria (gain and rotation) represent four 
initial parameters for adjustment: 


¢ assume the front current meter 1s correct and adjust the rear current meter in gain 
(RGAIN) or rotation (RTWIST), or 


¢ assume the rear current meter 1s correct and adjust the front current meter in gain 
(FGATN or rotator Fv Ts I), 
Because the sled instruments undergo significant stress in the surf zone, all sled transects 
conimenced with an orientation run outside the surf zone to allow for daily reorien- 
tation. A numerical orientation prior to each transect’s calculations was employed dur- 
ing data analysis. The resulting possible orientation adjustments for each transect and 
parameter (RGAIN,RTWIST, etc.) are shown in Table 6. 

The adjustments and resulting total radiation stress (S,/) were found to be linearly 
related, therefore an iterative method of linear interpolation. or “regula fals:’ method 
(Gerald, 1978) was utilized to determine the orientation required to minimize Sy be- 
tween the current meters. Convergence to within a radiation stress gradient tolerance 
of 10-3 Jm~* was generally obtained within six iterations. S,,(f) spectra with these ad- 
justments were perused to ensure that the strict numerical definition of a conserved ra- 
diation stress gradient achieved a realistic and comparative spectral shape for the current 
meters. 

The selection of orrentation by gain or rotation does have constraints when dealing 
with bi-directional radiation stress spectra, Which commonly occurred during SUPER- 
DUCK. Examples of typical bi-directional spectra are schematically illustrated tn 
Figure 15. The spectra indicate swell with a 10-second mean period from the southeast, 
While wind waves with a 5-second period and significantly more energy approached from 
the northeast. Spectra from the front current meter are assumed correct and are de- 


picted by a solid line, whereas the spectra from an incorrectly oriented or calibrated rear 
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current meter (and requiring gain or rotation to make its spectra coincide with the front 
current meter spectra) are depicted by a dashed line. In the first example, gain is the 
obvious correction Whereas in the second example, rotation is the obvious correction. 

All three current meters were calibrated in a flow tank before and after the exper- 
iment by different laboratories. The agreement in pre- and post-calibration was + 3% 
in gain and 0.01 — 0.04 ms7? in offset. This agreement is considered excellent for meters 
exposed to the robust surf zone. With this agreement, coupled with the subjective ana- 
lyses of the different spectral shapes corresponding to different orientation combina- 
tions, and the field accuracy of current meter alignment (+ 2”), it was decided to utilize 
Current meter rotation as the orientation criteria. Rear current meter rotation 
(RTWIST) was arbitrarily selected over front current meter rotation for the entire data 
analysis, except for Run 18-1 for which front gain was used. Resultant RTWIST angles 
are Within the values of field accuracy alignment. Sensitivity to this rotation and the 
effect of other rotations (FT WIST only or both FIWIST and RTWIST) is addressed in 
the error analysis discussion in Chapter VII. 

After completing this rotation for the rear current meter, the side current meter was 
rotated (STWIST) in a similar manner and its S,J values served as a first-order com- 
parison to the other two current meters values for S,/. These orientation rotations were 
performed after the @, coordinate transformation were fixed for the calculations inside 
the surf zone for that particular transect. This post-calibration procedure provided a 
means of improving the accuracy of the current meter calibration and alignment. The 
procedure 1s essential to obtaining realistic physical results from the data and 1s a dis- 


tinguishung feature of this research. 


FE. MEAN CURRENT OBSERVATIONS 
1. Overview 

Three-dimensional portraval of the mean current structure 1s illustrated in Fig- 
ure 16 and Appendix A (Figure 35 through Figure 40). All current vectors represent 
34.1-minute means. The numbers indicate the data run location. The short dotted lines 
above each run location indicate the mean water level (MWL) measured by the sled’s 
pressure sensor for that location and time. The different arrowheads delineate the height 
of measurement as illustrated in the three-dimensional box. MSL is indicated by the 
long dashed line. Only current meters which were completely submerged for the entire 
run duration are depicted. Thus some run locations have only one or two current vec- 
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Figure 15. Constraints on current meter rotation or gain correction for bi- 


directional spectra outside the surf zone. 


Table 6. INSTRUMENT ORIENTATION PARAMETERS. Rear and front cur- 
rent meter gains and rotations (“TWIST”) required to conserve the radi- 
ation stress gradient outside the surf zone. 
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Rear Gain(%) Front Gain(%) Rear Twist (deg) Front Twist (deg) 

=), 25 5.54 ee 0.80 
1.04 -1.03 -0,29 
rae | 4.38 -2.86 
Love! 0.58 0.12 
5720 -4,94 =i Jal 
1.40 -1.44 0.02 
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e ‘.’ indicates no iterative convergence or acceptable spectrum shape. 


e ‘Twist’s’, or current meter rotations, are in units of degrees. 





Wessletecas first deploved oi 1] October during the first of two SLPERDLCK 
storms. Hugh surf (/7,,, greater than 3 m) prevented sled deployment beyond the surf 
zone, therefore there was no instrument orientation run for that day. The high surf also 
prevented completion of a minigrid survey by the CRAB. Therefore the bathymetric 
profile, as well as the instrument orientation parameter was obtained from the 12 Octo- 
ber data. The unrealistic c, values later calculated for 11 October highhght the necessity 
for an instrument orientation prior to a transect. The 11 October data set was therefore 
deleted from further consideration. 

The 12 October data revealed a meandering rip current. Two transects were 
completed with Runs 12-1 and 12-4 being the orientation runs outside the surf zone. 
Run 12-4’s data were not used because of a Zeiss system malfunction. Therefore Run 
]2-] 1s the only data run outside the surf zone in Figure 36. The strong offshore flow, 
O(0.5 ms~!), at all three levels above the bed was substantiated by comparison with 
current vectors from the longshore array. Although this rip current eliminates homo- 
geneity for Sf calculations, it provides a unique data set for a rip current study separate 
from this research. 

October 15 1s the first dav with an established bar-trough svstem. Run 15-1 is 


the orientation run outside the surf zone. Strong flow, O(1.0 ms7!), is apparent inside 
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Figure 16. Mfean current structure across the bar on 16 October 1986. Numbers 


represent data runs, short dotted lines are MWL’s, and distance 
alongshore is arbitrary. 
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the surf zone with a 50% reduction in speed outside the surf zone. All speeds are nearly 
uniform in depth with slight speed reductions for sensor elevations closer to the bed. 
Multi-level current directions in the trough and outside the surf zone are co-linear. 
Angular deviation of currents with height (maxiunum 20°) are apparent at the location 
of maximum breaking, with the lower sensor elevations exhibiting a bias toward offshore 
flow near the bed. 

October 16’s data have the strongest velocities of all the sled data with speeds 
exceeding 1.0 ms~!. Run 16-2 was the orientation run outside the surf zone. Current 
directions and speeds were nearly uniform with depth, with current speeds exhibiting a 
slight decrease closer to the bed. Strongest flows were on the offshore side of the bar. 

Significantly weaker currents were measured on 17 October. Run 17-1 was the 
Orientation run outside the surf zone. The three trough runs all produced negligible 
current. A combination of low tide and low wave heights [ //,,,, O(0.2ms7!) J on this day 
may have caused most waves to break on the bar and thus cause the negligible current 
in the trough. The onshore flows at Run 17-4 and 17-5 were due to the current meters’ 
location in the upper water column of the shallow water. 

Strong flows, O(1.0 ms~!), were again measured on 18 October. Runs 18-1 and 
18-6 were orientation runs outside the surf zone for two separate transects. Run 18-6 
Was the start of a second transect which was abruptly terminated when the sled mast 
snapped in the high surf during Run 18-7. A nearly uniform honzontal velocity gradient 
(eh/ex) is observed with a maximum current located offshore and just before the bar 
crest. Vector uniformity is observed. 

Two transects were conducted on 2! October. Runs 21-1 and 21-6 were orien- 
tation runs outside the surf zone. Unfortunately current velocities during both transects 
Were too weak and variable in direction for meaningful analysis. Hence this data set was 
eliminated from further consideration, leaving only four (15-18 October) of the original 
nine transects for further analysis. 

2. Discussion 

Profiles of longshore currents across the barred beach generally showed signif- 
icant spatial variability. Maximum flow ts observed at the steepest offshore slope of the 
bar, generally just below the bar crest, while weaker currents are found outside the 
breaker zone and in the trough. A typical distribution of longshore currents across a 
barred beach is illustrated in Figure 16. Longshore currents are generated after wave 


breaking (Run 3 and 4) and are extended outside the surf zone, possibly by turbulent 


ay) 


momentum exchange. The vertical distribution of longshore and cross-shore velocities 
are compared to normalized water depth (Figure 17). Current measurements for one 
run are denoted by a continuous line connecting only those current meters which were 
entirely submerged for the entire run. Thus some shallow locations have onlv one or two 
current means illustrated. The plot of October 16 longshore velocities serves as a legend 
and has numbers beside the mean current symbols. The numbers represent the chrono- 
logical order of runs in a transect, with the most offshore run being labeled as number 
l. 

The magnitude of the cross-shore flow 1s generally much less than the longshore 
flow and is of O(0.2 ms~!). The longshore flow velocities on three of the four data days 
are of O(1.0 ms~'). The cross-shore velocities exhibit a relative onshore flow in the upper 
third of the water column and relative offshore flow in the lower third of the water col- 
umn. The lower mean currents after wave breaking tend to flow more offshore, thus 
exhibiting a limited “spiral” distribution with height and indicating the presence of an 
undertow flow (an offshore flow near the bed). The observations are in agreement with 
cross-shore theorv (e.g., Svendsen et al., 1987). The longshore velocities indicate nearly 
depth uniform flow with a slight velocity increase with elevation. This depth uniform 
flow confirms the appropriateness of the depth-integration assumption for longshore 
flow. Longshore vertical current structure did not exhibit a logarithmic profile, therefore 
the profile method for the determination of bed shear stress coefficients is not applicable. 

3. Statronarity and Homogeneity 

Stationarity over the transect was investigated by examination of S,/ calculated 
from measurements outside the surf zone by a fixed-bed pressure and velocity (puv ) 
sensor designated “south tripod”. The formulation for this calculation is discussed later 
in equations (78) and (79). Sequential 34.1-minute spectra and resulting Soe were Cal- 
culated for 15-18 October (Figure 18 and Appendix A (Figure 41 through Figure 43). 
The short lines between the arrowheads at the top of the Su graphs represent sled run 
durations, with the numbers above the lines identifying the sled run number of that day. 
These comparisons are considered a strict indication of stationarity because of the sen- 
sitivity of SF to changes both in wave height and wave incidence angle. 

The ‘puv’ data were collected by the FRF for four-hour periods, therefore these 
data were not completely synoptic with sled deployment times. The ‘puw data were 
digitized at 2 Hz. Spectral analysis was accomplished bv Bartlett ensemble averaging 


which resulted in 31 degrees of freedom. Chi-squared confidence limits for 95% confi- 
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Figure 17. Nlean cross-shore U and longshore velocity V profiles versus normalized 
water depth for 15-18 October. The 16 Oct V profile serves as a legend 
for the data symbols, where the numbers denote chronological order 


of runs in a transect with the most offshore run labeled as number 1. 
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dence yields a range for spectral estimates (S,,(/)) between 0.66 and 1.81 of the spectral 
estimate calculated. Ninety-five percent confidence limits for the integrated or total ra- 
diation stress (SJ, see equation SQ) were based on the Fisher-Z transformation (Miller 
and Ereund, (oss) 


Z=+ in| | (49) 


l—r 





where the correlation coefficient (r) is: 
——— (51) 


and the covariance 1s calculated by the integrating the co-spectrum of u,v. 

The October 15 runs inside the surf zone (15-2 through 15-4) appear to have a 
reasonably constant (within the bounds of the confidence limits) S,/ value with no more 
than a 10% variation about the mean. Both October 16 and 17 data sets have a constant 
Sf with minimal variability. October IS data also appears to have a reasonably con- 
stant S,7, with the exception of Run 18-5 which occurred at the time of an S,/ increase. 
Fortunately, this run is within the trough, where flow conditions are not as rapidly af- 
fected by offshore change as are locations on the offshore side of the bar. Even with this 
S,/ increase, the variation is no more than 20% about the mean. 


Tidal variation during an individual data run was approximatelv 5% of the 


y 
mean Water depth. Tidal variation over an entire transect ranged from 14% - 50% of 
the mean transect depth, with the longest transect duration having the largest tidal var- 
1ation. The tidal range for each transect is illustrated by the short dotted lines in 
Figure 16 and Appendix A (Figure 35 through Figure 40). 

Another test for stationarity 1s illustrated in Figure 19 where //,,, and wave in- 
cident angle («) calculated from the ‘puw’ sensor are plotted from 0500 - 2000 E.S.T. for 
each day and compared with sled run times. H,,, values are denoted by an “x” and a 
values are denoted by an “o”. Each short line at the top of the graphs represents the 
duration of a sled run, with the number above the line denoting the sled run number of 
that day. Comparison of Figure 19 with Figure 18 and Appendix A (Figure 41 through 
Figure 43) offers some insight into the nature of radiation stress. For example, the de- 
crease in S.J observed on 15 October from 1000 - 1200 was due to a slight reduction in 
wave height and a significant decrease in wave incident angle. In conclusion, no signif- 


icant change in S,/, H,,, , or o was observed during the times of sled deplovmiame 
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Figure 18. Offshore S,, spectra and S,/ from ‘puy’ gage 621 for 16 October 
1986. Numbers indicate data runs. Intervals are 95% confidence. 
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Stationarity requirements are therefore assumed satisfied. It should be noted that the 
Stationarity being investigated here 1s measured over relatively long temporal scales 
(measured in hours) whereas the temporal term examined later in Chapter VII is meas- 
ured over much shorter temporal scales (measured 1n minutes). In other words, the as- 
Sumption of stationarity at this stage of the data analvsis does not preclude 
contributions from the temporal term (e.g. due to edge waves) to the momentum bal- 
ance. 

Similarity of mean current vectors in the experiment area was considered a test 
for homogeneity inside the surf zone. Similarity was considered satisfied when all cur- 
rent directions were in the same quadrant with speeds + 50%> of an overall mean for the 
same relative on-offshore position. This definition was chosen considering the bounds 
of instrument accuracy, current Vertical variability, the number of instruments in the 
experiment area, and the variable elevations of the longshore array current meters. No 
visual evidence of rip currents was observed durmg 15-18 October. Two-dimensional 
plots of sled-deternuned current vectors were created for each transect, along with the 
fixed-bed array of mean current vectors, temporally centered for the sled transect. These 
comparisons were then overlayed on the coincident bathymetry (Appendix B). The 
bathymetry allows for isolation of cases where bathymetric forcing is dominant. Vari- 
able current direction is observed in the trough on 17 October. However this 1s not 
considered violation of homogeneity because of the extremely weak current speed in the 
trough on that day. Perusal of this data indicates the assumption of homogeneity for 
15-18 October is reasonable. Data sets from 15-18 October are considered to meet the 


dual assumptions of homogeneity and stationarity. 
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Figure 19. H,,, (‘x’) and wave incidence angle (‘0’) versus time for 15 - 18 October 
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1986. Numbers indicate data runs. 
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VI. SURFACE WIND STRESS (t#) 


A. INTRODUCTION 

The objective of this chapter is to determine representative surface wind stress for 
the surf zone based on atmospheric drag coefficients calculated from surface wind stress 
measured outside the surf zone and wind speeds measured by the sled anemometer inside 
the surf zone. Dr. S. SethuRaman and an investigative team from North Carolina State 
University calculated surface wind stress for the end of the FRF pier during SUPER- 
DUCk. Wind speeds Were measured with a three-component Gill anemometer at a 
height of 18.7 m above mean sea level (MISL) and wind stress was calculated using the 
eddy correlation method. Although this wind stress measurement was outside the surf 
Zone, it 1s a simultaneous and uncontaminated measurement over Water and within 5U0 
ni of the sled transect areas. 

A synopsis of the synoptic scale meteorological conditions during October 1986 is 
Irsted in Table 7. Sled deplovments were on 11-12, 15-18, and 21 October. Predonmiiem 
meteorological features were onshore winds and an unstable atmospheric boundary laver 
(as indicated by the air-sea temperature difference in Table 7). 

The Gill anemometer has three polystvrene propellers which are rigidly mountee 
with their axes at 90° from each other. These propellers drive three dic. generatam 
whose output 1s directly proportional to wind speed. The anemometer has a threshold 
of 0.1 ms—! . a distance constant of 0.8 m. and can be sampled at 100 Hz. Although 
additional wind measurements Were acquired at 22 m and 13.7 m, the eddy correlation 
method is selected over the profile method for wind stress deternunation due to its m- 
creased accuracy. Air temperature (7,) and relative humidity (RH) measurements were 
acquired by a YSI 44018 Thermuilinear thermistor and a slow response relative hunudityv 
sensor (Model 207) manufactured by Campbell Scientific, respectively. Both instruments 
were mounted at 18.7 m above MSL. A VSI temperature probe measured sea temper- 
ature (7,) at a depth of 0.6 m below MSL. 


Wind stress determined by the eddy correlation method is computed by: 
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Table 7. SUPERDUCK SYNOPTIC-SCALE METEROLOGICAL DATA (OC- 
TOBER 1986) (from SethuRaman et al., 1987), 


Date Wind Field Brief Description ipa louds Remarks 
l Soe 3-7. High dominates 6.4 Clr 
2 Nat. 3-7 High dominates a5 Cu.As S-B Hazy 
3 SE-SW 1-6 High dominates Lae Ac.As S-B 
+ SW 7-10 High donunates de Clr 
5 SW-NE 6-7 Frontal passage ore Cu,Sc S-B 
6 NW-NE gust 13. Post-frontal region -2.] CUA Sma sivell 
7 NE-SE 3-5 High dominates -3.5 Ac.Cu, B-S 
§ SESW High domunates 2 Sc.As B-O 
9 Evie-olre 1-5 Cold front approach Oe Cu,Sc S-B Rain 
10 hie 13-16 Cold front passage =e SieSc Oren 
1] Melt — 16 Post-frontal region -0.8 Suse O-b Sct 
he NE 6-10 High dominates -0.5 Ac.As § 
13 SE 3-8 Cyclogenesis along shore 1.0 Cu.As B- Rain 
14 SW 6-§ Cold front approach 1.4 iGrest-O Kain 
PD ee 7-i1 Frontal passage -6.8 Seu Swell 
16 NE-NW 7-13" Post-frontal region -3.3 Stise O-B Swell 
1/7 NE-NW 4-& High dominates -2.3 Cu.Ac § 
18 NW-NE 9— 14° High dominates -3.7 nrs-6 Swell 
19 NE-NW 6— 13" High dominates -2.9 Cu S-B 
20) NE 5-8 High dominates -2.0) Clr 
| Re-St 1-5 High dominates -1.3 Clr 
7a SW 3-5 High donunates he Clr 
a3 SW 4-6 High dominates URS, Cs $-B 
24 Be - NE 5-9 High dominates 0.0 Ac,Cs $-B 
5 NE-E 7-10 Warm front section 0.0 seic © iain 
26 Be-os) O-/ Cold front approach 4.1 St.CirO 
a] p= Sl: 3-3 Cold front passage gee Sc.As B-O Rain 
28 NE-SE 2-8 High donunates -0.6 Cu 
i) Sywei-5 Transition of 2 highs ney Sce.Cs B-O 
30 SW-NE 4-7 Front sweeps through ale Sc.Cs B-O 
31 pee pust 12 High dominates -0.§ Meee 5-1} 
Comments: 
ime scattered, B: broken, Clr: clear. O: overcast 
Bee clogenesis occurred in nearby region on IIth and 16th 
3. Strong upper-level trough accompanied cold front during 14-18 Oct 
4, * indicates high-energy cases of SUPERDUCK experiment 


S. 1, ~ 1, was measured at 1400 EST which was the approximate mid-time 
of the sled transects. 


where u.v,w denote wind speed in the horizontal and vertical direetions and the primes 
(‘) denote deviations from the mean. Hourly means of wind stress data (SethuRaman 


et al., 1987) allow determination of an atmospheric shear velocity (4, ): 





Us = ( a ie (52) 


Where t 1s the surface wind stress. 

Inherent in the w determination is the inclusion of atmospheric stability. The de- 
scription of the turbulent regime of a temperature-stratified friction layer is based on the 
Monin-Obukhov similaritv theory, which defines the stability-dependent vertical wind 


speed gradient as: 


atk Ls z a 
es oe (2) 
“eHene 2 = height above the surface (m), 

L = Monin-Obukhov length, 
K = von Karman’s constant (0.4), and 
Qn = a stability function. 


The non-dimensional ratio z,L 1s a measure of atmospheric stability with negative 
values denoting unstable conditions and positive values denoting stable conditions. z, L 


is defined as: 


Z Ske Ms 54) 
= (3 
L if iis 
Mere ao = virtual temperature (degrees Kelvin), and 
ee = virtual temperature scaling parameter. 


7, is often considered to be the mean virtual potential temperature within the surface 


laver (Geernaert et al., 19874) rather than its more strier deitmipionmea: 


T,=7,(1 + .6lg¢) (5 


where q = specific humidity. In this research, 7, in degrees Kelvin was substituted for 
ip 


° 
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Ivectatine equation (00) toi the Suriace (Zz = 2,, where u = 0) to an arbitrary 
height (Z = z,u = u,), an expression for a stability-dependent logarithnuc wind velocity 


profile is obtained: 





e. 


ee ( In = _ ral — )) (56) 


Determination of z’L requires an iterative process, since 7\s 1s also a function of 
stability. With all z L’s initiallv set to zero, and the humidity roughness scale (=,,) and 
teniperature roughness scale (z,,) set equal to 10~°, an iterative approach (Appendix C) 
based on the stability flux profile relation of Businger et al. (1971) is utilized to determine 


ziL. YY, and ‘Y, are assumed equal and formulated as (Large and Pond, 1981): 


ay foie 
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The air at the air-sea boundarv is assumed saturated. 
The stability function. @,,, 1s described by the Businger-Dyer formula for unstable 


Stratifications (Businger et al.. 1971) as: 
Pn =U mo) (38) 


Mievexpression for ‘Y,, is (Paulson,1970): 
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which is an empirical representation of comprehensive observations over smooth terrain. 


B. DETERMINATION OF AN ATMOSPHERIC DRAG COEFFICIENT FOR THE 
SURF ZONE 

Atinospheric shear velocity and L are independent of height in the surface boundarv 
Javer (Panofsky and Dutton, 1984) and are assumed to apply for the sled transect lo- 


cations in the surf zoue located within approximately 500 m of the surface wind stress 


(Sas 
tsa 


measurements. Atmospheric shear velocity and L for the sled runs are linearly interpo- 
lated from the hourly pier data for the temporal midpoints of each sled run. 

Boundary laver investigators universally convert drag coefficients to be appheable 
at 10 m height to allow intercomparison between different investigations. For equation 
(56) to be valid, a stabilitv-dependent wind speed at 10 m (Uj) is required. Ujy9 1s ob- 
tained from equation (36) as: 


. Lbs <m 10 =m 
— —_—_— =e —_—=— —_ — a i} 
c 10 Us K [ In( 10 ) 35 ara if ) Hen (60) 


Uy9 1s first calculated from the hourly pier data. Then two drag coefficients for a 
height of 10 m at the pier are determined -- the stabilitv-dependent Cy), (from equation 
I] with Uj) — U, substituted for W) and the stability-independent or neutral drag coel- 
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The surface boundary condition ( UW, ). represented by a mean current. was taken 
account to determine the wind speed relative to the moving water surface. C, 1S assumes 
zero for the end of the pier. U, 1s represented by the mean surf zone current for tiieeaes 
locations. 

Wind, wave height, atmospheric stability, and initial pier Cj data versus time for the 
seven davs of sled deplovment are shown in Figure 20 and Appendix A (Figure 44 
through Figure 46). Drag coefficients and pier and sled wind speeds are adjusted to 10 
m elevation. Sled anemometer winds are represented by the first. mid, and last two- 
minute mean Wind of the run. The clear arrowhead depicts the hourly mean wind vector 
for the pier wind and the solid arrowhead depicts a 34-minute mean wave direction vee- 
tor for a location near the pier end. The point of the arrowhead indicates the time for 
Which the vector 1s valid. Wind and wave directions are relative to the shoreline which 
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Figure 20. Wind speed, wave height, atmospheric stability, and initial Cy data for 
October 15-16, 1986. Beach is assumed parallel to the v-axis. Open 
arrowheads denote wind direction and closed arrowheads denote wave 


incident angle. Winds are adjusted to z = 10 m. 
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The absence of hourly 10’L (z L stabihty indicator for z = 10 m) data and wind data 
indicates no boundary laver information was available. Absence of hourly 10 L data 
when a comcident wind speed is depicted indicates a 10. L value more negative than the 
limit plotted. 

Boundarv layer stability ts a primary concern In interpreting the wind data. There 
were no stable atmospheric conditions during SUPERDUCK. C9 was plotted so as to 
remove the effect of stability in an effort to relate the drag coefficient variation touam 
other physical parameter, such as wave height, wind direction, or wind speed. Cayy9 18 
the standard drag parameter utilized for wind stress cormnparison between different in- 
Vestigators and serves as a basis for examining wave influence on wind drag. Neither 
wave height. wave direction, nor wind direction exert obvious influence on the drag co- 
efficient as shown in Figure 20 and Appendix A (Figure 44 through Figure 46). 

Increased instability places the atmosphere mm a convective instability or local free 
convection regime where buoyancy forees become equal or greater in maenitude than 
mechanical mixing forces. Wvyngaard and Cote (1971) conducted measurements of ver- 
tical velocity and temperature under unstable conditions and observed local free con- 
vection to occur at z L < -0.5. Using the Businger et al. (1971) stability Mux promie 
relation, Garrett (1977) showed that as instability increases, not only does C, increase; 
but roughness becomes significantly more important. The Businger ct al. (1971) iterative 
approach 1s mapplicable in a local free convection regime. Therefore an empirical ap- 
proach is deemed more appropriate for surf zone C, determination under localeine 
convection. 

Sled runs of z L < -0.5 generally occurred during periods of low onshore wimdisjpeces 
Therefore the drag cocfficrents in Figure 20 and Appendix A (Figure 44 through 
Figure 46) depicted for low wind speeds and very unstable conditions are inaccurate. 
Onshore winds of low speed and diabatic (unstable) conditions are representative of a 
sea breeze. A sea breeze is a local coastal wind caused by the temperature difference 
when the sea is colder than the adjacent land. It usually occurs on relatively calm and 
sunny davs (Huschke, 1959). Review of surface weather analysis maps (Appendix D) 
and wind speed. wind direction, and stability (Table 7 and Table 9) indicate sea breeze 
conditions for most of the sled data run periods. 

Measurements of surface wind stress above the swash zone under the effect of a sea 
breeze were mide by simultaneous wind and temperature profiles over a shore transect 


near Fort Walton Beach. FL (Hsu, 1970). The sea breeze in the surface boundary sieen 
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was determined to be in the atmospheric free-convection regime. A least-squares fit to 


Me data (igure 21) vields: 
Cry X 10° = 10.06 —1.30(Uyy — U,) (63) 


where Uj) is the mean wind speed at 10 m above the sea surface. It is noted that the 
drag coefficient under the sea breeze condition 1s higher than that under more stable 


conditions which is in agreement with e.g.. Roll (1965). 
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Figure 21. An empirically-derived stability-dependent drag coefficient obtained un- 


der a sea-breeze condition (from I]su, 1970). 
To determine appropriate drag coefficients for the sled transects based on the 
SLPERDUCK wind stress data, the following procedure 1s emploved: 


Smee) te Sled s anemomieter fuiled (Runs 18-1, 18-2, and al! runs on 21 October), 
a linearly-interpolated hourly mean Uj, from the pier is used as the sled Ujo. 


e A hinearly-interpolated Cz, from the pier is assigned for the temporal mid-point 
of the sled run. 


e When the actual z,L was more negative than -] or the wind speed was less than 
6.5 ms~!, C9 1s determined from the empirical relation of equation (63) where the 
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sled’s Uj) 18 determined from equation (69) and the z represents the run-varving 
height of the sled anemometer above niean water level for that particular run. 


e¢ When z L was more negative than -] and the wind speed greater than 6:5 ms —™ 
Ciao 1s linearly interpolated between the value preceding and following the one in 
question. This is felt to be a better approximation than not using stability. 
C. ADJUSTMENT OF C, FOR THE SURF ZONE ENVIRONMENT 
Since the surface wind stress was measured outside the surf zone, a C, adjustment 
is required to more realistically reflect the surf zone environment. Early investi@augme 
of oceanic C,’s sought a dependence on wind speed. Subsequent investigators (e.g., Hsu. 
1974, 1986a; Byrne, 1982: Geernaert et al.. 1986) began to search for a sea Statemeee 
pendence for C,, since surface gravity waves are assumed to act as roughness elements 
(z,) for the momentum transfer process. Most of these investigations were done for open 
ocean conditions and therefore are not directly applicable to the surf zone. 
1. Uniqueness of surf zone air-sea interaction 
The uniqueness of surf zone air-sea Interaciion as Opposed to open ocean air-sea 
interaction can be summarized bv three factors -- shoaling of the bottom. wave refrac- 
tion, and wave breaking. 
a. Shoaling effect 
\fomentum transfer is based on the relative wind speed difference between 
the absolute wind and wave speeds. During shoaling. there is a wave steepening and 
celeritv decrease (relative wind increase) of surface gravity waves as thew sisi 
Geernaert et al. (1987b) state that C, increases with decreasing values of the ratiovaaas 
where c, 1s the phase speed of the domimant locally-generated gravity wave. Thus, larger 
than expeeted drag coefficients are observed over the sea during shoaltie COilmdinieme 
(Smith, 1930: Geemaentct al. 19se. 
b. Refraction 
In most cases, the wave propagation direction 1s at some angle (up to 
1$0°) to the predominant wind direetion due to refraction of the waves mn shoaling water. 
Over the open ocean, C, 1s found to be larger if the wind ts blowing at an angle different 
than the wave direction (Large and Pond, 1981; Smith. 1980). Geeraert et aieGie aa 
proposed wind stress enhancement in advance of a cold front due to surface gravity wave 
convergence. Swell generated behind the cold front would travel through the front and 
interact with Waves generated ahead of the front and propagating parallel to the front. 


This interaction would result in a chaotic sea of larger wave amplitudes, steeper gravity 
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waves, and more wave breaking -- very similar to a robust surf zone. This enhancement 
was observed by Davidson et al. (1988a, 1988b) during the FASINEX 1986 experiment. 
c. Wave Breaking 

Moderate to high surf conditions undeniably result in increased surface 
roughness, which would result in increased drag coefficients. Although some work 
(Amorocho and DeVries, 1980; Melville, 1977) has quantified relations between breaking 
waves and increased drag coefficients, these investigations have dealt with the open 
ocean where the onset of breaking waves is produced by a dramatic increase in wind 
speed. 

2. Surf Zone C, Adjustment Formulation and Results 

The wind stress data from the end of the FRF pier is assumed to include both 
the shoaling and refraction effects on C, as one transitions from the open ocean to the 
Bemezone, Ihe mcrcased C, effect due to shoaling from the end of pier to the surf zone, 
as Well as the additional refraction O(15 — 25° maximum). is considered to be small 
compared to the breaking wave effect. 

The approximate mean distance from the swash line to the end of the pier 1s 500 
fieeieach dav at 0700 EST, the distance from the mean swash line to the outer extent 
of the surf zone was “stepped off along the pier. These rough values for surf zone width 
on sled deployment days are listed in Table 8. The surf zone never extended to the end 
Siete pier, although on I! October, the farthest surf zone extent nearly did. With 
onshore winds, the increased roughness would not take effect until one was a 
Maeasuieawle distance within the surf zone. Therefore, the increased rougliness due to 
the breaking waves of the suf zone is not reflected in the end of the pier wind stress 
data. 

The purpose of this research is not to explicitly characterize the C, for the surf 
zone, but to intelligently select an appropriate C; with which to address longshore sur- 
face wind stress as measured by a single anemometer in the surf zone. The increased 
surface roughness due to surf zone breaking waves may be conservatively likened to the 
chaotic seas observed bv Davidson et al. (1988a, 1988b) in advance of storm fronts. 
From a composite of wind stress measurements taken from multiple chaotic sea cases 
found in advance of storm fronts, the drag coeflicient was found to be increased by ap- 
proximately 33% over that expected for non-chaotic seas and equal wind speed 
(Davidson et al., 1988a, 1988b). Therefore an increase of 33°% in the drag coefficient 


(determined from equation 61) was deemed appropriate for all runs except those runs 
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Table 8 SUPERDUCK SURF ZONE WIDTH FOR SLED DEPLOYMENT 
DAYS (from Crowson et aly 1ve0 


October Surf zone width (m) 
1] S00 
}2 320 
i) 120 


16 100 
17 SO 
18 SO 
| 100 





occurring under strong diabatic conditions where equation (63) was uulized. Hsu's 
(1970) relationship is empirical and therefore the C,; values determined by equation (63) 
would alreadv have the increased roughness considered. Drag coefficient values for each 
run and procedural step are listed in Table 9. AH the values are of O(107%). 

Surface wind stress is determined from equation (11) using the C,s hsted in 
Table 9. Atmospheric density is assumed constant at 1.25 kgm7? and @ is obtaimed from 
the end-of-prer wind direction. Uj) — U, 1s substituted for W where Uj) is the stability- 
dependent mean windspeed measured by the sled anemometer for the entire data run 
(approximately 34 nunutes), corrected to an elevation of 10 m. U, 1s the mean surf zone 
current measured by the sled current meters. The term U,)—U, 1s obtameadmirom 


vectorial subtractron. 
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VIE. LONGSHORE MOMENTUM BALANCE 


The tongshore momentum balance is calculated based on field measurements of 
wind, waves, and currents acquired across a barred beach. Spatiallv-dependent bed 
Shear stress coefficients and the relative rmportance of the various terms in the momen- 
tum balance are calculated. Wind force to wave force ratios are calculated from the field 
data and compared to the results of the wind forcing model in Chapter III. This chapter 
concludes with an error analysis of these calculations. 

The stope arrav did not perform as required for determination of the wave-induced 
longshore momentum flux (Sy) in equation (1). Pre- and post-cahbration values for the 
Slope array differed as much as 10%. Assuming that breaking wave turbulence outside 
the surf zone is negligible, the gain required to match the slope array’s radiation stress 
(Syx(A) to the current micter’s radiation stress (S,,(/)) was O(30 - 60%). Therefore, slope 
array data were judged in error and were not used in this analvsis. Studies addressing 
turbulent momentum flux were suspended. 

The slope array loss necessitates a rewriting of the longshore momentum balance 


equation as: 








a VX = Eee He “7 (64) 


The turbulent momentum flux gradient 1s combined with the wave-rnduced momentum 
flux gradient to form a combined momentum flux gradient as measured by the current 
meters. The S,/ term implies integration over frequency (equation 80). Determination 
of surface wind stress is addressed in the Chapter VI. Calculation of the remaining terms 


In equation (64) will be addressed separately in this chapter. 


A. TEMPORAL, OR ACCELERATION TERM (6AZ,/ér) 
The time-averaged depth-integrated mean momentum per unit area mn the longshore 


direction (.\/,) 1s composed of a steady (M,) and unsteady (.!/,’) flow component: 
YF ? >) ell ” ? lad 
Wy = We Vi pliner ip pv’ az (65) 
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The steady component is calculated every 64 seconds with the term, 7 + A, being taken 
as the mean water level (MWL) measured by the pressure sensor for that period of time. 
This resulted in thirty-two mean values of 64-second duration for A, . 

The unsteady flow contribution to the mean momentum, .\/,’, is attributed to the 


mass flux of the wave motion: 


n 
M, = |, ov dz (66) 
Where v’ is only Known at z = QO from wave theory. By expanding in a Taylor series 
about z = 0 and retaining only the first order term, 47,’ may be expressed independent 


of height in the crest-trough region: 


" 
i = |, pv'(0)dz = pv(O)n (67) 
Using linear theory and substituting for v’(O) and »: 


aa (68) 


pe l l 
My = > Psa 


where a is wave amplitude. Rewriting equation (68) in terms ofa finite frequency band, 


a significant wave incident angle (a(f)) , total wave energy (E(f)), and wave celerity (C(f) 








yields: 

q? fr E(f) ae. 

M,/ = I, Cif) sin a(f) df (69) 

where 

a= —— tanh kh | (70) 
w(f) 


Throughout this analysis, wave number (k) is calculated non-iterativelv by the formu- 
lation of Wu and Thornton (1986). 

Spectra of “significant wave angle” (Appendix F) are calculated as (Higgins et al., 
1981): 


1s 


Su 
[EV nN] (71) 


sin a(/) cos o.(f) = 


The energy spectrum is calculated from velocity measurements. Time-averaged 


depth-integrated kinetic energy (KE) may be written as: 


Ae Re a eee 
KE =0.5 p _ Au + vay jaz (72) 


where w 1s the component of the flow velocity in the vertical direction. Jo the finsmgmemes 
of approximation, there 1s equipartition of energy between potential energy and kinetic 
energy for linear progressive waves (Kinsman, 1965). The total energy, E(f), is calculated 


as twice the kinetic energy: 
E() = 214,()1° [6,10 + GN] (73) 
where 


e( sinh(2kh)] 


| Hy (f) | = 2 
Sk[ cosh"k(h + z,)) 


(74) 
and G, and G, are the one-sided spectral density functions for the horizontal velocity 
components in the x and y directions and z,, is the measurement depth. Calculations 
of equations (69) through (73) indicated that the time-varving momentum of the waves 


for this data set was one-to-two orders of magnitude less than the time-varying mo- 


mentum of the steady flow. Therefore the temporal term may be simplified as: 








car AVin, 

ar ~ AI 2) 
where 
Wie Nh = pV an) (76) 
and AA = time step over which the finite difference 1s calculated, and 


4 + h = MWL as sensed by the sled s pressure semsmrr 


A value of hdl is first determined for each ensemble of 64 seconds in the data run. 


Since the temporal scale of the longshore current is unknown, a range of ensembles 
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(comprised of 1 to 11 segments of 64-second duration) are defined using a running mean 
interval, Ar’. Then a finite difference is calculated froin these larger ensembies for a 


specified Ar. Finally, a mean temporal term 1s calculated: 


k 
aM, ok BG a 
Ar &k Be |e ) 
PS 


for the k increments of A‘/,/Ar. This mean temporal term is assigned as the temporal 








term for the sled run. Calculations of the mean temporal term (equation 77) with and 
Without a running mean, as Well as variation of the time step (with subsequent change 
in the k increment), made negligible difference for this data. Therefore Ar = 64 seconds 


is used. The selection of Av’ is discussed in paragraph D below. 


B. RADIATION STRESS GRADIENT (0S, 7 / 6x) 
Longshore momentum flux, as defined in equation (2), 1s calculated from current 


meter measurements and applving linear wave theory: 

2 
En) 7 | H(/)| Cn (78) 
where | H,(/)|2 is a linear transfer function: 


oL sinh(2kh) + 2kh] 


| H(A\° = 
a 4kL cosh*Ath + z,,)] 


(79) 


and C,,, is the real part of the cross-spectrum of u’ and v’. It is assumed in the applica- 
tion of the transfer function (equation 79) to the measured covariance that the vertical 
distribution of turbulence in the water column is the same as in the waves, and is de- 
scribed by the transfer function which is a function of depth and frequency 
Meteure 22). 

The total radiation stress (S,/), where total indicates frequency summation, is de- 


temmuned from: 
Sy 
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Figure 22. Linear wave theory transfer function used to calculate S,,(f) from cur- 


rent velocity. showing decreasing frequency dependence with depth 


across a barred beach. 
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where f, and f, represent the lower and upper frequency limits of the sea-swell band un- 
der investigation and Af represents the resolution of the spectral estimates. S,/ is the 
quantity required for the momentum balance of equation (64). 

Analysis of the longshore component of momentum flux (S,,(/}) was accomplished 
using Fast Fourier Transforms and ensemble averaging. Ensemble averaging requires 
segmentation of a time series into several shorter time series of equal length, and calcu- 
lation of a sample spectrum for each of these sub-series. The mean of these spectra 1s 
a smoothed spectral estimate. The advantages of ensemble averaging are minimal com- 
puter storage and computational ease. 

The data were initially demeaned and detrended. The 34.1 minute data runs were 
segmented into 32 sub-series of length 64 seconds. The S,,{f) used in equation (80) is 
the Bartlett smoothed spectral estimate (Jenkins and Watts, 1968) from 32 sample 
seectia, All spectral estunates (S,,, 7, etc.) have 63 (v = 2m - 1) degrees of freedoin, a 
Spectral resolution (Af) of 0.015 Hz, and Chi-squared 95% confidence limits equal to 
0.73 and 1.45 times the spectral estimate itself (Jenkins and Watts, 1968). Since sea and 
swell are of interest, lower and upper frequency hmits are established at 0.06 and 0.44 
Hz, respectively. These limits are used whenever frequency summation was required for 
anv parameter, 

Radiation stress spectra may be negative at one frequency while positive at another 
frequency, as in the case of sea and swell approaching from different quadrants. S.,(/) 
spectra transitioning from positive, through Zero, to negative, and vice versa, may not 
be accomodated with a log-linear plot, therefore S,.(/) spectra are plotted on linear plots. 
Related spectra are also plotted in this manner for consistency. 

The calculation of the gradient of radiation stress requires numerical differencing 


over an appropriate length scale: 


(51) 


The spatial step (Ax) was initially designated as the distance between the front and rear 
current meters after adjustment for sled yaw angle. A radiation stress gradient (equation 
81) determined with this definition of Ax is labeled to have been determined by the “lo- 
cal” method. The Ax distance on the sled was approximately 2 m and the “local” finite 


differencing scheme is: 


ds 


Ey: T E 

OSyx (Soe Jrear ci. (Spe Vicon clm 

——— 2 ee (82) 
Cx os 

A second approach to the determination of the radiation stress gradient, called the “in- 

tegrated” method, utilized the distance berween sled run positions for Ax. This distance 


was O(20-40 m) and this finite differencing scheme Is: 


T T. ig 
Le - (Sy Ja ft (Sie )y = (Syx 3 
Ox 4 Axi ogee 





yx yx yx 


(‘So _& y= 
Axy_4 AXx,_3 


(83) 


where the numbers | through 4 denote the two pair (front and rear) of current meter 
data available from two sled run positions (see Figure 23). The estimate of the S.J 
gradient is based on four independent current meter locations covering slightly different 
transect segments and slightly different Ax’s. These segment and Ax differences are 
generally less than 10% of the distance between sled mast locations. It is therefore as- 
sumed that the four estimates represent four sample means of the true S,/ gradient be- 
tween the two sled locations. From the Central Limit Theorem, the standard deviation 


of sample means 1s written as: 





Tose 
G Ox 
+7, = Or CASS ak aa (84) 
: mdi _—_— ett 
\ ae N 


Where n is the number of sample means. Thus the standard deviation of the sample 
mean is reduced by 50% with the use of four sample means rather than one. This fact 
will be used in determining the measurement error for 6S,, ‘/ Ox in paragraph E. Se- 
lection of the “local” method versus the “integrated” method proved to be of significant 


importance. 


C. SEA SURFACE ELEVATION (7) 
Although not specifically required for the momentum balance, sea surface elevation 
spectra provided a useful reference against which to compare radiation stress spectra. 
Sea surface elevation spectra (G,(f) ) are determined from linear wave theory using 


data from both the sled’s pressure sensor and current meters: 
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Finite differencing scheme for 6S,, '/¢x. from “integrated” method. 


2 P 
GLY) == | 17;()) | Gp) (85) 
where 


CO sh?(kh) 


| 2,()|° = ——————— 
: cosh*Lk(h + z,,) ] 


(86) 


and G, (f) is the one-sided energy density spectrum for the pressure head (p,) measured 
by the pressure sensor. The pressure head 1s determined from the hydrostatic equation: 


P — Pam 
Ph = pg (87) 


where p is the absolute pressure sensed by the pressure gage and p,,,, 1s the atmospheric 
pressure: 


The formulation for determining G,(/) from current meter velocities is: 


GN = |HAAI CG, + 6 ($8) 
where 

gz S] 2k 
HAL? = sinh(2kh) ($9) 


2ek cosh’ Lk(h +z) 


G,(f) produced by these two formulations are in good agreement for the three 
verticallv-stacked current meters and the pressure sensor, thereby providing confidence 
for the 7 spectra. Sea surface spectra produced in Appendix F are formulated from 
equation ($5). 


Wave height parameters are calculated from pressure-derived 4 spectra: 
Mea 
ims = (8 g,,)2 (90) 


| (Gil ) 


where the total variance 1s expressed as: 


f 
=) Gl af (92) 


Jy 


D. RESULTS OF THE MOMENTUM BALANCE 
Field measurements of wind, waves, and currents are used to calculate the local 
longshore momentum balance at various locations across a barred beach. [From the 
momentum balance, bed shear stress coefficients and the relative importance of the 
various terms comprising the balance are calculated. Mean values are determined rela- 
tive to the bar. Wind force to wave force ratios are calculated from the field measure- 
ments and compared to the results of the wind forcing model in Chapter ITI. 
I. Bed Shear Stress Coefficient 
A spatially-dependent bed shear stress coefficient (c,) may be determined from 
equations (18) and (64) as: 
OS 6M 
ot tL 
SS (93) 
pluly 





> 


Resultant spatially-dependent bed shear stress coefficients are listed in Table 10 and 1l- 
lustrated in Figure 24. The column titled “c{/)" are the cys determined utilizing the “lo- 
cal” method for the radiation stress gradient. These coefficients are representative for 
the sled run position itself. The surface wind stress and temporal variability terms were 
34.1-nunute means taken only from the sled run being investigated in anticipation of 
gaining insight into the temporal variability of all terms in the longshore momentum 
balance. A running mean of 12.8 minutes was selected for determination of both the 
temporal term and the radiation stress based on radiation stress statistics of Freilich and 
Pawka (1987). They showed that, given n independent realizations of a time series of 
u(t) and v(t), that the Gaussian approximation for true S,, statistics is valid for n > J] 
independent realizations. An enseinble length of 64 seconds and a minimum number of 
ensembles (12) was selected for the analysis. Unfortunately, this procedure resulted in 
significant cf/) variability. Utilization of a 34.1-minute mean 6S,, 7 / 6x made negligible 
difference in resultant c{/). Two possible reasons for the large c(/) variability may be 
that the radiation stress 1s much more variable across the surf zone on a small spatial 
scale than previously imagined and the Ax is too small to resolve the difference between 
two large numbers (see equation (82) and paragraph 2 of the Error Analysis section in 
this chapter). 

The columin titled “c{7)" are those cys determined using the “integrated” method 


for the radiation stress gradient. These coefficients are representative for the portion of 
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the bathymetric profile between sled run positions. The wind stress and temporal vari- 
ability terms are the 34.1-minute means from the two appropriate sled run means. The 
“beach” position along the transect is defined as the intersection of the bathymetric 
profile with the MWL as determined by the sled pressure sensor on the run closest to 
the beach. Extremely weak longshore currents inshore of run 17-6’s location explain the 
variable cfi)’s. They are disregarded. The bed shear stress coefficients determined by the 
“integrated” method exhibit less variability than those from the “local” method and are 
all of O(10-3). The “integrated” method is considered to be a reasonable method for the 
determination of the radiation stress gradient. 
2. Relative Contributions of the Longshore Momentum Balance Terms 

Relative contributions of wave forcing (¢S,//éx), wind forcing (77), and 
temporal variability of mean momentum (6M, / €t) to the total momentum balance are 
determined by first defining their total contribution (TOT) as the cumulative sum of 


their absolute values: 


a ip 
VX 


Gs 
ONG ai a! | 
COX ‘ 





(94) 








Their individual relative contributions (rc) are then determined by dividing by the total 














contribution: 
Sy. 
Cx 
wave, = TOT (95) 
| Ea 
wind,, = TOT (96) 
Cur 
OM, Cr 
ar Je aoa et) 


Relative contributions from the wave, wind, and temporal terms, as formulated 
in equations (95) to (97) and utilizing the “integrated” method, are calculated. Results 
are presented in order of consecutive runs (Table 11) and bv relative location along the 
transect with respect to the bar (Table 12). Relative locations are delineated as “off- 


shore” of the bar, “on the bar” and immediately before the bar, and in the “trough” re- 
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Table 10. SUPERDUCK BED SHEAR STRESS COEFFICIENTS (OCTOBER 
15-18, 1986) 


October 15, 1986 October 17, 1986 
Run Run CZ) 


17-1 0.002 
17-2 0.003 
17-3 0.004 


17-4 0.026 


17-5 0.099 
17-6 0.119 
October 16. 1986 October 18, 1986 


pat ee 


Run c/{2) Run C2) 


16-2 0.001 
0.003 

16-3 
0.004 

16-4 18-3 
0.00] 

16-5 18-4 
0.001 

16-6 18-5 
0.002 

Beach Beach 


Reeindicates c-vaiue < 5. x 10 
cf) are determined by the “local” method 


c(i) are determined by the “integrated” method 


relative contributions determined by the “integrated” method 





gion. The mean wind stress remained nearly constant throughout each transect of the 
surf zone. However, the temporal term, as averaged over 34.1 minutes, was variable in 
sign and magnitude. The location of maximum wave breaking along the transect can 
be identified (Table 11) where wave forcing’s relative contribution is maximum. Two 


peaks in wave forcing correspond to two breaking areas - waves breaking offshore of the 


bar and waves reforming and then breaking on the nearshore beach face. The relative 
contribution of wind stress inside the surf zone generally increases in value as one pro- 
ceeds onshore from the location of maximum breaking, due to the decrease in wave 
forcing. The relative contribution of wind stress can then decrease again inshore at the 
location where waves re-break on the nearshore. The relative contribution of the 
temporal term appears to be largest in the trough where wave forcing is reduced and 
possibly where the effect of longer period waves (e.g., surf beat or edge waves) is more 
pronounced. 

The mean relative contributions for the wave, wind, and temporal terms (0.86, 
0.08, and 0.06, respectively) are essentially identical for both “offshore” and “on the bar” 
(Table 12). However in the “trough”, the mean relative contributions (0.64, 0.17, and 
0.19, respectively) are significantly different. The relative importance of wave forcing 1s 
reduced by 25%, whereas the relative importance of surface wind stress is doubled, and 
the temporal term relative importance 1s tripled. 

A comparison is made between the wind force to wave force ratios produced by 
the wind forcing model in Chapter 11] to those determined for the SUPERDUCK data. 
The wind to wave force ratios in Table 3 and Table 4 are based on a planar beach, 
therefore only the offshore side to the SUPERDUCK bar (tanf ~ 0.03) is consideneas 
The SUPERDUCK force ratios are obtained from ratios of the relative contributions 
of surface wind stress to radiation stress gradient in Table 12 for the “offshore” and “on 
the bar” categories. The wind force to wave force ratio for 15-18 October fielayda is 
0.09. The mean wind speed was 9.6 ms~' and the mean wave incident angle for the 
“offshore” region was 12°. Using linear interpolation in Table 3 and Table 4 results in 
a wind to wave force ratio of 0.12. Thus the field measurements are in good agreement 
with the results of the wind forcing model in Chapter I1I, with differences attributable 
to non-fully arisen seas and the presence of swell froma different direction than the wind 
waves. 

Bed shear stress coefficients are calculated for cases without the surface wind 
stress and without both surface wind stress and the temporal term. The percentage 
change in cy ranges from near Zero to 100%, with a mean change of 10-30%. Thus under 
certain conditions for this data set, the surface wind stress and temporal terms can be 
of first order in the momentum balance. 

A mean cy based on relative postion along the transect with respect to the bar 


1s determined. Bed shear stress coefficients determined from the “integrated” method are 
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Table 11. RELATIVE CONTRIBUTIONS OF THE LONGSHORE MOMED- 
TUM BALANCE TERMS FOR SUPERDUCK 


Relative contributions determined by “integrated” method. Run number indi- 
cates offshore run of the required run pair. 


Runs 17-6 through 17-8 had negligible currents of variable direction. Their cy 
results are disregarded. 
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Table 12, RELATIVE IMPORTANCE OF MO- 
MENTUM BALANCE TERMS FOR LO- 
CATIONS RELATIVE TO THE BAR 






















Offshore 


On the Bar 


Trough 


ist 0 
0.77 0.17 0.06 
























assigned and a mean cy determined for each of these relative locations (Table 13). The 
results indicate a decrease in mean c, as one proceeds from offshore to onshore with 


mean cy values equal to 0.004 (“offshore”), 0.002 ("on the bar”), and 0.001 (trough’). 
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The shoreward decreasing cy was initially thought to be the result of decreasing bed 
roughness. As previously stated, bed roughness is a function of both sand characteristics 
and bed form. Beaches generally exhibit a shoreward gradation of sand grain size from 
large to small. Waves tend to sort the sand with the largest sand residing in the most 
energetic area of the surf zone. This gradation would be commensurate with the maxi- 
mum cy found just inside the surf zone and a decreasing cy further onshore. Unfortu- 
nately, field measurements of bed roughness were overlooked during the SUPERDUCK 
experiment and cannot be used to substantiate this hypothesis. Interviews with 
SereROUCK divers were conducted several months after the experiment, however 
their information about bed roughness was not specific enough to draw conclusions. 
An October 1982 sampling of sand grain size across the FRF beach showed shoreward 
increasing grain size (Birkemeier et al., 1985). It 1s reasonable to assume that grain size 
gradation would not change over a few years, therefore bed roughness is discarded as the 
reason for the c, trend. Field and laboratory studies have shown that bed shear stress is 
enhanced when the near-bottom flow field is composed of mean current and oscillatory 
feemwave) Components (Grant and Madsen, 1979, 1982; Grant et al., 1984). Thus 


shoreward decreasing wave action would result in shoreward decreasing Cy. 


Table 13. MEAN BED SHEAR COEFFICIENTS DETERMINED FOR LO- 
CATIONS RELATIVE TO THE BAR 


rods!) On the Bar Offshore 


15-1] 
16-3 
16-2 
17-1 
ler = 2 
18-1] 
18-2 


CD GO ~~) ON GN CA Ch 
i 
ip Jibs (On mom Wen ow) JL. 


l 
| 
l 
] 
l 
l 
l 


Mean Cr 





The variation in cy due to varying Av’ in the temporal term) was examined by 


analyzing sample groups where Av’ was specified as a running mean of (approximately) 
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I, 5, 10, 12, 14, 18, and 22 mintites. The mean vananonss cy due to varying Av’ for the 
SUPERDUCK data was meciictole. 


E. ERROR ANALYSIS 
An error analysis 1s conducted to assign appropriate error bounds on the resultant 
spatially-dependent cy values. Approximated errors are the use of a linear wave theory 
transfer function to integrate the velocity covariance measurement over the water col- 
umn and errors in stationarity due to tidal variation. Three sources of error which are 
specifically quantified are incident wave angle error, propagation of measurement error, 
and finite differencing error for the radiation stress gradient. 
1. Incident Wave Angle Error 
Radiation stress (equation 3) is sensitive to the incident wave angle (a). The 
determination of an accurate « 1s dependent on several other angles. The Zeiss elec- 
tronic surveying system Was required to measure the sled vaw angle. The Zeiss error 
(Table 14) is assumed negligible in the « determination. The adjustments (Chapter V) 
to transiate the sled évstrtnent’s three-dimensional coordinates into the FRF coordinate 
system, and the measurements by the portable Zeiss svstem which made them possible, 
are considered to have negligible error. Three angular sources of error which cannot be 
neglected are beach angle orientation error, current meter orientation error, and current 
meter directional gain error. 
a. Beach Angle Orientation Error 
An analysis is conducted to assess sensitivity to beach angle orientation er- 
ror. The beach angle orientation 1s calculated from a mean bathvmetric contour deter- 
mined in a least squares sense as described in Chapter V. For this analysis. 
monochromauc, linear, unidirectional wind waves are assumed for a straight-and- 
parallel planar beach. Using Snell’s Law, the wave incident angle (05) at the inshore 
current meter may be determined given the wave incident angle (@,) at the offshore cur- 
Fent Micter: 


_j| sina, C; (98) 


where wave celerity (C) is determined from equation (70). The radiation stress gradient 
(from equation 3) in the cross-shore direction may be stated as a finite difference be- 


tween radiation stress calculated at the two on-offshore current meters Where the sub- 
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Figure 24.  cyand percent breaking wave (Q) values as a function of offshore distance 
(October 13-18, 1986). (cy is above and Q is below the profile line). 


vl 


script “rear” refers to the rear or offshore current meter and front to the irommsen 


inshore current meter of the sled: 
CSyx —  LEnsin(a + Ao,) costa + Ao.) rear — LE nm sin(a + Aa.) cos(a + Aa.) Irom (99) 


Cx Ax 





where Ax, is the beach angle orientation error. Total energy (E) is calculatedirom: 


| 
E= wedi (100) 


Root-mean-square wave height (/7,,,,) for monochromatic shallow water waves can be 
related to depth (h) in the inner surf zone by equation (39). For SUCPERDUCKiiamas 
of /7,.., to h ranged from O(0.3) outside the surf zone to O(0.4) inside the surf zone. 
The radiation stress gradient (equation 99) versus orientation error for rep- 
resentative depth, slope, and incident wave angles 1s illustrated in Figure 25. The three 
angles plotted for each depth are centered about a tvpical incident wave angle encount- 
ered during SUPERDUCK for that particular depth. For example, an incident wave 
angle of 25°, a depth of 3.5 m, and a beach orientation angle error of —5°, would result 
in an overestimation of the radiation stress gradient by 10%. On the other hand, a ra- 
diation stress gradient measured in a shallower depth with a reduced incident wave angle 
(h = 1.5m, a = 15°) and the same error angle of —5° would result in an overestimation 
by 30%. Note that the overall governing parameter 1s incident wave angle. Therefore. 
for a fixed beach orientation angle error, the relative radiation Stress error increases with 
decreasing wave angle, and hence with decreasing depth. For SUPERDUCK, errors in 
mean beach orientation angle are estimated to be less than + 2° and 1.0 mi Wage 
minimum depth in which data were acquired. Therefore radiation stress gradient errors 
for 3.5, 2.5, and 1.5 m depths would be + 5%, + 7-8%, and + 10%, respectively ies 
radiation stress error can be equated to approximately a 10% cy error, as long as wind 
stress and temporal variability of mean momentum are small. As the magnitude of wind 
Stress and temporal variability increase, the effect of a fixed radiation stress error on ¢¢ 
becomes less. Varving beach slope from 0.01 to 0.10 or varying Ax from 2 m (local” 
method) to 25 m (“integrated” method) had negligible effect on the relative radiation 
stress gradient error due to an incorrect beach orientation angle. Therefore cy, error due 


to inaccurate beach orientation angle 1s estimated to be less than 10%. 
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Figure 25. Sensitivity of CS,, / ¢x to beach angle orientation error for varying inci- 


dent wave angles and depths and tan 8 = 0.03. 
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b. Current Meter Orientation Error 

Orientation of the current meters axes to the sled frame axes was considered 
to be within I — 2° (Chapter V). Assuming conservation of the radiation stress gradient 
outside the surf zone, all current meters were individually oriented by a previously de- 
scribed selection process (Table 6). This numerical rotation 1s considered to have elim- 
inated any CUprent meter onentationemor 

To assess the effect of other rotations which minimized the radiation stress 
gradient and which had realistic spectra shapes, rotation subsets within Table 6 combi- 


nations were also examined. For example, 


e inrun 16-2, rotating both current meters by a lesser amount instead of just the rear 
current meter by —2.86°, or 


¢ in run 15-1. rotating the front current meter by 0.27° instead of the rear current 
meter by —0.29°, 
yielded no substantial difference in the magnitude of c, along that transect. In other 
words, by using the concept of conservation of S,, to initially orient the current meters, 
the resultant c,s Were insensitive to the manner in which the orientation was done. 
c. Current Meter Directional Gain Error 
Current meter directional gain error is the error introduced by using a single 
gain and bias factor to calibrate current flows for all directions. The current meters were 
initially calibrated in a laboratory flow tank where gains and biases were determined for 
meter probe angles of 0°, 90°, 180°, and 270°. These angles correspond to the four 
sensor probes protruding from the current meter ball. Using the gain and bias factors 
for the u and v sensors determined for these four probe angles, the current meter was 
then checked against a controlled flow of O(0.5 ms!) for meter probe angle increments 
of 15° . The standard deviations from the controlled flow for sled current metersuieees 
and 4 were 0.012, 0.012, and 0.009 ms7!, respectively. This error is included within the 
measurement error for current velocities addressed in the next section. 
2. MNfeasurement Error 
To ascertain the effect of direct measurement error on the final determination 
of cy, an error propagation analysis is conducted for both the “local” and “integrated” 
methods for cy determination. Standard deviations (a) of directly and indirectly meas- 
ured variables used in equation (93) were estimated or determined from laboratory call- 
bration results (Table 14). 
The coefficient of variation (Barrv, 1978), or fractional standard deviation 
(Young, 1962), is defined as: 
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Table 14. MEASUREMENT ACCURACY 


Direct Measurement 
parameter 
water densitv 
Zeiss vert control! 
Zeiss horiz control! 
c/m velocity 


aS 


> 


tvpical value 
1022 kgm~? 
reaon 

40 m 

u = 0.25 ms! 
v = |. ms7! 


Be lela 


time 

Wind stress + 0.0- Colina 
MWL (Paro) 0. 2m 

wind speed 5 (ons! 
atmos. press. 1010.0 mb 
atmos. temp. ‘ 


Indirect \feasurement 


atmos. drag coefl. 
atmos. density 


e where + refers to an approximate standard deviation. 


' from Birkemeier and Mason (1984) 








V.== (101) 


where V is the arithmetic mean of replicate measurements of a variable. For the coef- 
ficient of variation determination (Table 14), the strict Y has been replaced by typical 
(order of magnitude) values, which are listed. 

Calibration of the sled current meters was performed before and after the ex- 
periment by two different calibration laboratories. Pre-calibration gains and biases are 
listed in Table 15. Post-calibration gains varied by + 3.0% and 0.01 - 0.04 ms~! in bias 
or offset. A conservative estimate of measurement standard deviation for the current 
velocities 1s + 5% to include the unmeasured environmental source of error between field 
and laboratory conditions as well as the uncertainty (+ 3%) associated with using a 


single gain factor for all frequencies (Guza and Thornton, 1980) and directions. 


The general rules of measurement error propagation can be stated for four 


arithematic operators as (Barry, 1978): 


ay 2 2 2 
Osum or diff = mY, C4 T (Cpe 9 Cn ae na (102) 
CA \2 CR .2 Lae? i 32 5 
Goroduct = x ABC ... n/(“4) + (=) ~~ (=) + .. + a (103) 
Gquotien = th /(tyP + (BY (104) 
quotient = Re \ A B 
SO power = +nA”'c, (105) 


where the standard deviations a, og etc., refer to the + error of the parameters (e.g., 
A. B) beins analyzed (young. 12 a2): 


Table 15. SLED CURRENT METER CALIBRATION 


Cit Sled Sled Pest Correlation 
Ser,= Letn dist Gai Channel) Vel Gane Cocfl (mg 
9118 rear, mid l 6 us 0.951) 0.000 0.99944 
5 vy 0.916 -0.008 0.99970 
S762. front’lower J i WW" (lSoe Ors 0.99899 
a es Oe 0.99878 
$354. side upper D S us 1.457 -0.030 0.99949 
. Vv e457) SOU, 0.99955 


¢ current meters were Marsh-McBirney, Model 551 
¢ calibrations by U.S.G.S, Gulf Coast Hvdroscience Center, NS Wie 


e formulation: velocity in ms—' = (gain x voltage) + bias 


The wind stress measurement error has the largest V, of all the parameters in 
Table 14. V. can reach as high as 40%, depending on the manner and conditions under 
which the measurement was taken. Conditions during SUPERDUCK were conducive 
to accurate measurements of wind stress. These conditions were: 
¢ over-water measurements taken during onshore winds, 


e atthe end of a 500 m pier, 
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e ata height of 18.7 m MSL which 1s 11.1 m above the pier deck, 
e using the eddy correlation method, and 


e using a Campbell CA27T Sonic Anemometer as backup for the Gill anemometer. 


Standard deviation of atmospheric density (p,) 18 estimated from equations 
(103), (104), and the equation of state for a perfect gas (Haltiner and Martin, 1957). 
Therefore, using a subjectively-determined and conservative value of + 30% for the wind 
Stress measurement error, a C, coefficient of variation of 30% is determined from 
equations (11) and (103). 


Rewriting equation (93) as: 





. a 
ONee . OM, 
Be TT “A 
= ee ee (106) 
plu,ly plu,ls plia|v 
or 


the measurement error in cy, or the standard deviation of Ch is the result of the additive 
os of the three terms on the right-hand-side of equation (106), abbreviated as A, B, C 
in equation (107). Applving equations (102) to (104), results in standard deviations for 
the four terms in equation (107) as listed in Table 16. 

Standard deviation determination of the radiation stress gradient (¢,,), requires 


o of (OS, 2 / 6x), which in turn requires o of S,7. SJ or S,, are calculated from equation 


) 
(78) using the directly measured parameters u’,v’,p, andh. An upper bound for o of 
Se can be formulated by approximating for shallow water waves: 


S) 
Se = ph |, Cm = = Col, = Pit 0G, (108) 


where coy,, 1s the covariance of the u and v velocities and is bounded by (e.g. Davis, 
1986): 


es Oy Oy (109) 


Therefore o, 7 can conservatively be approximated as: 
yx 
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ogra ST [Ce (Sy (Ee + OY (110) 


og, and a, are negligible compared to o, anda, (Table 14) and may be neglected. 


One may then write: 


(111) 


(112) 





The standard deviation of Ax is much less than the standard deviation of AS,J and 
therefore may be neglected. The standard deviation determined by equation (112) for 
the “mtegrated” method 1s reduced by a factor of one-half to account for that method's 
increased sample size (see equation 84). Standard deviation of the radiation stress gra- 
dient for each of the sled’s data runs, using both the “local” and “mtegrated” methods 
for cy determination, were calculated. The mean for radiation stress gradrent of all four 
transects is listed in Table 16. 


The standard deviatron for the wind stress (t?) was determined from equations 


(11) and (103), and Table I4. Standard deviation of o\f, / or was determined by as- 





sumung: 
om, P AM, — AlpDV) (115 
cr ee Aes Al 


since VW < ie and then applying equations (102) through (104) and Table 14 values. 
A typical mean cy value of 0.003 was used in equation (101) to determine the cy 
coefficient of variation m Table 16. Because of the relatively large incident wave angles 
and the measurement precision utilized in the SUPERDUCK experiment, propagation 
of measurement error resulted in reasonable cy values for both methods (Table 16). The 
“integrated” method's O-, Was superior to the “local” method because of the larger rela- 
trve error for the radiation stress gradient (equation 112) experrenced when Ax is small. 


The mtegrated method's error bounds for the cy range in Table 13 are + 4-15%. 
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Table 16. ERROR PROPAGATION ANALYSIS 


itemmifegn 107) o_ il cy Method 


A 0.0015 local 
0.00013 integrated 

B 0.00003 

C 0.00007 


Cr 0.0015 : local 
0.00015 integrated 


where tern) 


oaks 


sll 





3. Finite Differencing Error 
Thornton and Guza (1981) examined the error involved in representing an ar- 
bitrary continuous function by discrete measurement points. Using a Tavlor series ex- 


pansion for S,,: 





M6) hen GSne. 
a yx\ Ae YX 

eye + Ax) = Sy(x) + Ax ca “1 7 ee sen (1 14) 
and then solving for the true radiation stress gradient: 

ay os ado a3 
eee A Ax Os (115) 
Cx AX 2! Ax? 3! leg i: . 
or 

CS aS 

yx VX 
—. -;, 7 Error (116) 


where the sum of the second and infinite terms on the right-hand side of equation (115) 
is the error of the finite difference space approximation. 


Relative error (R) mav be defined as: 


oe 


Ree Error (117) 


AE 
OS) 


Cx 





Solving for R where only the first error term in equation (115) is considered, S,, for 


monochromatic shallow water linear waves inside the surf zone may be expressed as: 


S,, = E sina cosa x AD (118) 
Where 
ieee NE 
le ¢ PS) De 
sale oF ee 
D = Ape h, aiid 


> 
I 


complicated constant 


Substituting equation (118) into (117): 


Ay tan fi (119) 


Sa : h 


Wheres 71, ( Phormtom and Guzasos 

This approximation (equation 119) is adequate providing Ax and tan f are rea- 
sonably small. However, Ax and tan # are not always small in this apphcation. There- 
fore R is redefined to include a sufficient number of terms on the right-hand side of 
equation (115) to allow convergence to within 5% of R. Resultant computations re- 
vealed that equation (119) was sufficiently accurate (within 5%) for mOStGIaaaae 
SUPERDUCK data. However, using the “integrated” method between the most inshore 
data run and the beach required additional iterations (up to five) to stabilize R. 

The relative error in equation (117) 1s biased in one direction because the term 


of largest magnitude is written as: 








e Se A 
Ax yx Cie 
ZO ay 
SS ee (120) 
Ones Ch 
Cx Ox 
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and assuming negligible bottom curvature. Therefore the sign of R 1s fixed as the same 
sign as beach slope. The finite approximation for the radiation stress gradient may be 
Damitten a's: 


AS, aSy. 
A = (+R) 








“~ 
OX 


where R represents the overestimation of 0S,,/ 6x for positive slopes and the underes- 
timation for negative slopes. Finite difference relative error for the four transects of this 
research are listed in Table 17 and Table 18. They are listed for both methods of cy 
deternunation. The larger error for the “integrated” method value between the most in- 
shore data run and the beach was due to the combined effect of increased Ax, decreased 
depth, and steeper slope, a// which combine to increase R. The relative error is signif- 
icantly reduced when using the “local” method due to its much smaller Av. “Integrated” 
method’s R ranges from -27% on the nearshore side of the bar to nearly 200% on the 
nearshore slope. However, except for the gross error on the nearshore slope. “inte- 
grated” method's relative error is less than 50%. This would cause c, error due to finite 
differencing to be less than 50% (assuming negligible wind stress and steady state con- 
ditions in equation 93), 
4. Combined Error 

The various sources of error are now combined into a total, or combined, error 
relative to the bar. To establish a conservative combined error, finite differencing error 
is assigned a range of error rather than a bias and unquantified errors are estimated at 
= 20%. Mean finite difference error for the trough does not include the gross error, 
O(200%%), between the most inshore sled run position and the beach. Therefore the final 
cy Value for the trough does not reflect the nearshore slope. 

The maximum combined error is determined by adding the four maximum pos- 
sible errors listed in Table 19. However it is not reasonable to define maximum error 
in this manner since each of the four errors have a range of uncertainty associated with 
them and they are uncorrelated. It is very unlikely that the errors would all be a maxi- 
mum and in the same direction at the same time (Barry, 1978). Hence, assuming the 
errors are independent of one another, the combined error is obtained from equation 102 
and results in final momentum balance cy Values of 0.004 + 0.0010, for offshore the bar, 
0.002 + 0.0006 for on top and immediately before the bar, and 0.001 + 0.0003 for the 


trough region. 
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Table 17. FINITE DIFFERENCE RELATIVE ERRORS 


October 15 1986 Kelative Enron 
Sled Run oe) tan fs MWLi(m) R(%) 


ks 
DS, 


2. 
oan 


lee? 
21g 


log 
24.24 


Peeks 
44.65 


be, 
42.03 


SS S99 99 S98 
OO Sa SSS 


Beach 


October 16 1986 Relative Error (R) 
16-2 GOOlss 


16-6 


Beach 
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Table 18. FINITE DIFFERENCE RELATIVE ERRORS 


October 17 1986 Relative Error (R) 
Sled Run Ax(m) tan £ MW L (m) @) 


17-1 1.99 0.0182 Bro E: 
34.60 0.0224 De 2 12. 
17-2 1.99 0.0397 JL Se pa 
26.80 0.0364 ZO) a besiey 
les, 1.99 0.0309 1.43 Bs 
9.54 0.0265 1.3] ihe, 
17-4 1.99 O02 21 ele 3 
4.09 0.0232 1.14 6. 
EES oo 0.0166 1.09 2. 
SESE -0.0153 125 -27. 
17-6 1.99 -0.0036 1.40 0. 
4.88 0.0109 1.48 S 
17-7 1.99 Oy AOD Ve 
6.27 0.0276 1.49 Sk 
17-8 ok) 0.0495 1.43 Ss 
20.06 0.07106 Ont eae 
Beach 
October 18 1986 Relative Error (R) 
18-] 1.99 0.0113 3.80 Q, 
48.33 0.0259 3.09 a 
18-2 ee, 0.0405 2246 Re 
18.85 0.0369 I-94 yas 
18-3 ieee 0.0334 1.49 oy 
21.60 0.0258 Nes 34. 
18-4 1.99 0.0083 1.13 l 
Diol -0.0096 P30 -15. 
18-5 1.99 -0.0223 1.48 -2. 
33.7] 0.0438 0.74 1906. 
Beach 
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Table 19. COMBINED ERROR RELATIVE TO THE BAR ON THE BED 
SHEAR STRESS COEFFICIENT AS DETERMINED FROM A MO- 
MENTUM BALANCE 


eee “a 


| beach angle orientation error —_| angle orientation error 






finite differencing error 
mean R 
Ave 


= R 
Ay error = I/(1 + R) 









final cy 0.001 + 0.0003 | 0.002 + 0.0006 | 0.004 + 0.0010 
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VUl. LONGSHORE CURRENT MODEL 


A. INTRODUCTION 

A first-order longshore current model for a barred beach is proposed. The model is 
based on balancing the wave-induced radiation stress gradient with longshore bed shear 
stress, turbulent momentum exchange, and longshore surface wind stress. Stationary 
wave conditions and straight and parallel bottom contours are assumed. The waves are 
described as random, but narrow-banded in frequency and direction. Model results are 
compared to field observations and bed shear stress coefficients are determined by 


model-fitting. 


B. MODEL THEORY 
There are no field-verified longshore current models for a barred beach. Therefore 
a model which has been field-verified for a planar beach was selected for this analysis. 
The model (Thornton and Guza, 1986; hereafter referred to as TGS6) is based on a 
simple balance between the gradient of the wave-induced radiation stress, the bed shear 
stress, and the gradient of turbulence-induced radiation stress. The model’s assumptions 
are stated in the previous paragraph. The longshore current model is based on a random 
wave height transformation model (Thornton and Guza, 1983; hereafter referred to as 
TGS3). Modifications to the TG&6 model are made based on the availability of addi- 
tional data from SUPERDUCK and with the purpose of improving and broadening the 
model's applicability to barred beaches. These modifications are the inclusion of surface 
wind stress in the momentum balance, retention of incident wave angle dependency, and 
the use of a breaking wave height distribution based on a slope-dependent y and an im- 
proved empirical expression for the breaking wave percentage as a function of H, h, and 
y. Whereas TG86 found turbulent momentum exchange to be unimportant for their 
planar beach topography, turbulent momentum exchange is shown to be required in 
order to model current Velocities in the trough region of a barred beach. The TG86 
longshore current model is briefly summarized below with modifications described later. 

I. Wave Height Transformation 

This portion of the model is an extension of work by Battjes and Jannsen (1978). 


The wave height transformation is conceptualized as an energy flux balance given by: 
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= 12, (i 2a 


womisis 1G cos is the shoreward component of wave group velocity. % is the mean inci- 
dent wave angle of the narrow-banded waves, and < ¢,> is the ensemble-averaged dis- 
Sipatron due to wave breaking. Analysis of SUPERDUCK data, as well as resulis from 
TGS83, indicated that the ensemble-averaged dissipation due to bottom friction < é> is 
neghgible compared with < ¢, >, and therefore < €-> 1s omitted m equation (122). 
Based on analvsis of field data from Torrey Pines Beach, CA, the Ravleigh wave 
height probability density function (pdf). p(/7), was shown to qualitatively describe the 


measured wave heights across that planar beach. The Rayleigh pdf 1s described by: 
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where H 1s wave herght. The Ravteigh distribution describes all waves. 

Onlv the breaking waves are dissipative, resulting in changes in the momentum 
flux, 1e. Wave foreing. Therefore a distribution describing the breaking waves aione 1s 
required. Wave herght probabrhty densitv functions Were calculated from field meas- 
urements taken at Ft. Ord’s Soldiers Beach. Monterey, CA. These pdis weresiseqmae 
more reahstically model random wave height transformation across the surf Zone. 

Within the surf zone, breaking wave height distribution Lp,(H)J is approxi- 


mated by wercghting the Ravleigh pdf: 
plH) = WH) pA) (124) 


Where the weighting function W(H) < I. It is noted the breaking wave distribution ts 
not a pdf. but that the area of p,(/7) gives the percentage of the waves that are breaking. 
A simple form of W(H) which approximates the fraction of waves that break from off- 


shore to Saturation conditions 1s: 
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for saturated wave breaking conditions. Based on the It. Ord data, a modification to 


equation (125) which more heavily weights the larger waves 1s given by: 


Eli 2 5 
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Wave breaking dissipation is modeled as a simple periodic bore: 
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where B is a breaker coefficient of O(1) and f is the mean wave frequency which is the 
peak frequency for narrow-band wave spectra. The B coefficient accounts for differences 
in various breaker types and is a function of the proportion of foam region on the 
@eeaxer face (Figure 26). The peak frequency in the SUPERDUCK 4, spectra is used 
for f. 





Figure 26. Periodic bore used to describe spilling breakers. (from Thornton and 
Guzael9s3). 


The ensemble-averaged dissipation due to wave breaking is calculated by 
multiplying the dissipation for a single broken wave of height H by the probaoility of 


wave breaking at each height Lp,(H) J, and integrating over all H: 
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3 
<f,> = eer Hp.) dH (1233 
By substituting from equations (124) and (127) and evaluating the definite integral, 
equation (129) mav be written as: 


= te J 
22 Pe ———————— (130) 
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Energy flux (EC, from equation 122), 1s approximated using linear theory for C, at fre- 


quency /. and energy 1s calculated based on the Ravleigh pdf: 


| er ae I 2 : 
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Substituting equations (130) and (131) into (122) vields a first-order differential equation 
for the linear wave herght transformation model: 
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A spatially-dependent wave height CH,,,,(x)] can then be determined for an arbitrary 
bottom profile across the surf zone by numerical integration. given an initial value of 


fl 


ms the H,,,, profile 1s used as input for the lonsshore eunpent calemiauen 


2. Longshore Current 
The longshore current distribution across the surf zone 1s based on conservation 


of momentum. Rewriting equation (3), wave-induced radiation stress may be written: 





= - sing ; 
Se = eC, COS (133) 
where 
: . NN 
Sins SIN %,; 
= constant = ——— (134) 
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by Snell's Law of linear wave refraction. The subscript “o” refers to any arbitrary lo- 


cation. Therefore the wave-induced radiation stress gradient may be written as 
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where < ¢, > is substituted from equation (122). Assumung steady state conditions and 


negligible wind stress, the longshore momentum balance 1s written as: 


. ra 
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TGS86 used equation (136) and model-fitting to solve for the horizontal profile 
of longshore current across a planar beach. They obtained a hierarchy of solutions by 
choosing combinations of linear or non-linear bed shear stress, with and without turbu- 
lent momentum exchange. Turbulent momentum exchange was found not to be im- 
portant for their planar beach topography. They concluded that a non-linear bed shear 
Stress without turbulent momentum exchange was their preferred choice for gencral ap- 
plications, based on its fit to their data. A mean Cr for their four davs of data was 0.006 
+ 0.0007, 


C. MODIFICATION OF MODEL BASED ON FIELD DATA 

Modifications to the TGS6 longshore current model are the inclusion of surface 
wind stress, incident wave angle dependency, and use of a breaking wave weighting 
function based on a slope-dependent » and an improved empirical expression for break- 
ing Wave percentage as a function of H, h, and y. The use of turbulent momentum ex- 
change 1s found to be important for modeling longshore eurrent velocities in the trough. 
The inclusion of surface wind stress and incident wave angle dependency are covered in 
the Model Results section. The other modifications are addressed next. 

l. A Slope Dependent ; 

Since the TG8&3 model was developed for planar beaches, y was assigned a con- 
stant value of 0.42 (equation 39) and is the ratio H,,,./h for saturated wave breaking 
conditions. However. Sallenger and Holman (1985) found y to be a function of beach 
slope based on data from several field experiments. They determined an empirical re- 


lation for a slope-dependent y based on their data: 
y = 3.2tan B + 0.30 (137) 


i eva 2 isecOmMpared to beac Slope (Fisure 27) for both the SUPER- 


DUCK data and the data used by Sallenger and Holman (1985). The solid line in 
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Figure 27 is Sallenger and Holman’s (1985) empirical relation for 7. Most of the 
SUPERDUCK data has the same magnitude of scatter as previous data sets about the 
solid line, with even negative slopes (found on the inshore side of the bar) in general 
agreement. The exception is for SUPERDUCK data sets located in the trough and on 
the foreshore slope (Runs 11-3, 15-4, 17-6, 17-7, 17-8, 21-4). They have a significantly 
lower value of H,,,./h and are indicative of non-saturated breaking wave conditions. 
Obviously a slope-dependent y relationship 1s required for a barred beach. Based on 
Figure 27 and the general fit of the SUPERDUCK data, equation (137) is Wsedeigmaene 
SUPERDUCK analysis. 
2. Breaking Wave Weighting Function 

The weighting function, W(H), used in TG83 was based on only four data runs 
for a single transect of a beach with a tidal plateau. To improve on the formulation of 
W(H), breaking wave pdfs are calculated for the SUPERDUCK data. Time series of 
pressure head. [p,(r)], are linearly transformed into (2) time series by first calculating 
the complex Fourier spectra of p,(z), then multiplying by a linear transform function 
LH fh]. band pass filtering between 0.03 - 1.00 Hz, and finally applying an inverse 
Fourier transform to the #({f) function. 

The transfer function [(f) (from equation 38) is held constant at the 0.5 Hz 
value for the interval 0.5 I1z to 1.0 Hz so as to include the additional higher frequency 
Valance [On the da. 


increasing Hf) values. Filtering 1s accomplished by zeroing out the Fourier amplitude 


calculations, but not to overly bias them with the exponentially 


coefficients of the undesired frequencies before applying the inverse Fourier transform. 
The enure 34.1 min record is transformed at one time. As an example, the initial three 
minutes of Run 15-2’s p,(1) time series is compared with its transformed y{/) series in 
Figure 28. The transfer function H(f) and the G,{/) spectra determined from the entire 
y(t) series are also included. Wave periods and heights are determined from the 7(7) time 
series by the down-crossing method. Root-mean-square wave height is determined sta- 


tistically from the calculated (1) time series as: 


ie § a o He |" (138) 


i=] 


The statistical H,,,. from the y(t) series are compared with the energy-bascauyaees 
(equation 90) determined from both the G,(/} spectra discussed in Chapter VII and the 


G,(/) spectra calculated from the (1) record (Figure 28). Root-mean-square wave 
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Figure 27. Ratio #7,,,./h versus beach slope for various data sets.  Curcled data 


represent trough and foreshore slope locations. 
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height from all three sources compared favorably Lwithin (10%)], with the difference 
indicating a measure of non-linearity in the surface elevation record. 

As discussed in Chapter IV, visual observations of wave dissipation (i.e., wave 
breaking) at the sled were electronically marked on the data set records. These marks 
are then superimposed on the calculated 7(t) time series to identify which waves were 
breaking during each 34.1 min data run. From this procedure, wave heights and periods 
for all waves, and a subset for just breaking waves, are obtained. Pdfs of all waves and 
breaking waves for 15-17 October are illustrated in Figure 29 and Appendix A 
(Figure 47 through Figure 49). Since the data were measured sequentially and not si- 
multaneously, wave heights are normalized by the mean offshore wave height to account 
for any variability in the incident wave conditions. The pdf's are normalized to the mean 
offshore wave height measured bv a WAVERIDER buoy located 6 km offshore. 
Equipment malfunction precluded breaking wave identification on other days of the ex- 
periment. 

Breaking wave data from SUPERDUCK, as well as the four data runs from Ft. 
Ord (TG83), are compared in Figure 30. The ordinate is the fraction of waves which 
are breaking (Q) and the abcissa is the ratio, //,,.,/yh, where y 1s slope-dependent and 
calculated from equation (137). Thus for saturated conditions, both the ordinate and 
abcissa are equal to unity. The breaking wave fraction is obtained by dividing the 
numiber of waves identified as breaking by the total number of waves. The total number 
of waves in a data run Was approximately 400. The Ft. Ord data are indicated bv an ‘x’, 
SUPERDUCK data acquired in the trough are solid circles, and all other SLUPERDIRIAE 
data are indicated by open cirlcles. Three breaking wave weighting functions are illus- 
trated. The dotted line is the TGS$3 weighting function (equation 127). The chain-dot 
line is from a theoretical expression (Battjes and Janssen, 1978, Battyes and Stive, 1985) 
where the cumulative probability distribution of all wave heights (breaking or non- 
breaking) 1s the Ravleigh distribution, but abruptly terminated at the breaking wave 


height: 
Hims 
g= | fing +1 (139) 


Neither of these two functions provide a reasonable fit to the data. However, a function 


which substantially improves the fit to the data is: 
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Figure 28. Example p,(t) and (1) time series with corresponding G,(f) spectra. 
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Figure 29. 
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Wave height and breaking wave height distributions (hatched area) nor- 


malized to offshore 77,,,, for 16 October. The Rayleigh distribution is 


the solid line and empirical breaking wave distributrons given by 


equations (127) and (140) are the dotted and dashed lines, respectively. 
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H, 
Oe wey = | 1+ tanh s| See — 099) 1 — exo} -[F (1.40) 


and is delineated by a dashed line. This function (equation 140) has no physical justi- 
fication but 1s merely an empirical expression derived to better fit the non-trough data. 
It is noted that none of the three functions fit both the trough and non-trough data. 
The right half of the right-hand side of equation (140) is taken from equation (127) and 
is used to more heavily weight the larger waves towards breaking. Observations show 
that the largest waves at a particular location are more likely to break (TG8$3). Thus 
skewing of the pdf is substantiated by the breaking wave observations (hatched areas) 
seen in the pdf figures. Very little wave breaking occurs until H,,,,/yh reaches 0.8, then 
the percentage of breaking waves increases rapidlv. The addition of the SUPERDUCK 
data suggests a more rapid increase in the breaking wave percentage within a smaller 
H»,;/yh band -- almost to the point of describing a “step” function indicative of all waves 
breaking at a fixed breaker point. 

The ratio //,,,,/yh in Figure 30 has a certain amount of error associated with It. 
The error in y is the result of an imperfect fit to the data in Figure 27. Tidal fluctuations 
@iso aiiect the # term. Thus there are runs which have /7,,,,./yh ratios greater than one, 
vet were unsaturated. 

It is also noted that during field identification of the breaking waves, breaking 
Wave percentages were visually estimated to be often in the range between 60-90%. Yet 
subsequent data analvsis reveals the percentages to be much lower. This may perhaps 
be simular to the phenomena associated with visually estimating wave heights from ships, 
where it has long been recognized that the human eve perceives an average of the highest 
1,3 wave heights and not the true mean wave height. This lead to the designation of ship 
wave height observations as “significant wave heights” and not mean wave heights. In 
a similar manner, SUPERDUCK visual observations of breaking wave percentages were 
overestimates. Thus visual observations of saturated wave conditions from other ex- 


periments may need to be reviewed. 


D. MODEL RESULTS 
1. Wave Height Transformation 
As input to the wave transformation model, H,,,, and f measured during the 
orientation runs outside the surf zone (Runs 15-1, 16-2, 17-1, 18-1) are used as initial 


conditions. Mean wave frequency is assigned as the frequency of the spectral peak in 
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Breaking wave fraction (Q) versus the ratio 7,,,,/yA. 


Figure 30. 
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sea surface elevation. Mean incident wave angle is then assigned as the angle for f from 
the significant wave spectra for that day. [G83 and TGS6 assumed small wave incident 
angles and set the cos @ term in equation (122) to unity. Because of the large incident 
wave angles during SUPERDUCKk. incident wave angle dependency is retained 
throughout the analysis by inclusion of the cos x term. A forward stepping scheme is 


Meed to evaluate equation (122): 
feercos 2), = (EC, cosa), + (<e,>),Ax con 


All variables in equation (141) are Known except B, which is determined by model-fitting 
the calculated H,,,, to the field-measured H,,,,. Model-fitting 1s accomplished by mini- 
muzing the least square error between observed and modeled data. Results of the 
model-fitting are shown in Figure 31 and Appendix A (Figure 50 through Figure 52). 
ier 7,,. 1s represented by an “x for each of the sled runs and is determined from 
pressure sensor measurements. The spatial step size, Ax, was varied from 0.5 m to 10.0 
Miewith nechcible effect on the H,,,. profile. Therefore Ax = 35.0 m was arbitrarily se- 
lected to allow suitable cross-shore definition along with reduced computational time. 

A constant B coefficient of O(1) is determined by a least square fit for all four 
days (Table 20). Field and modeled H,,,, show good agreement except for the trough 
area on 15 and 16 October. Root-mean-square wave height 1s undermodeled by about 
35.0 for these exceptions. Spatially-variable B coefficients were also deternuned which 
Bese xact fits to the field duta. However their drastic spatial changes in B appeared 
unrealistic since the breaker type was observed to be similar for most of the transect. 
except on the very nearshore slope inshore of the last sled runs. A constant B cocthiwient 
for a transect was deemed appropriate. 

2. Longshore Current Model Results 
Longshore current 1s modeled by including surface wind stress in the momentum 


balance equation of equation (136): 


sin &, ; 


- 


b 4 é : 
—f > = —-% + ty —- SS 
3 : Y 4 Cy 


° (142) 


and using various combinations of the following approaches: 
¢ linear or non-linear formulation of t°, 
e with and without turbulent mixing, 


¢ spatially-variable or constant cy and N. 
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Figure 31. Nfodel-generated H,,,, and V compared with field observations versus 


distance for invariant cy and B and including turbulent momentum ex- 
change, 16 October. Model results are solid lines and field observa- 


tions are “x”’s. EjJevation is referenced to MSL. 
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In all cases, the spatially-invariant B coefficients listed in Table 20 are used. As with the 
B coefficient determination, cy and N are model-fitted by a least square error criteria. 

A spatially-invariant cp with a non-linear (and presumably more reahstic) i and 
without turbulent mixing results in poor agreement between modeled and field-measured 
longshore current velocities (Figure 32 and Appendix A, Figure 53). Although the 
general shape of the observed longshore current profile is modeled, there are large dis- 
crepancies between observed and modeled current velocities, especially in the trough re- 
gions. Current velocities in the trough region are in error by as much as a factor of 7. 
Order of magnitude for c is similar to the momentum balance determination of cy 
(Chapter VI1) for all 4 transects. 

Allowing cy to vary spatially provides some improvement in the model-fit. 
However major discrepancies still occur on the nearshore slope where a steep increase 
in velocity 1s modeled along with sigmificant underestimation of velocity in the trough. 
The spatially-variant cy solution is dismussed because of these discrepancies. 

Although TG86 deternuned turbulent momentum exchange to be unimportant 
for their planar beach data, turbulent momentum exchange 1s investigated for this anal- 
ysis as a possible mechanism to transfer momentum across the bar and into the more 
placid trough region where wave breaking is reduced. Turbulent momentum exchange 
is expressed as the last term in equation (142) and 1s approximated using an eddy 


viscosity formulation: 





— Lp AV 
Syx = . puvdz = — péD a (143) 
where D 1s total depth and the eddy viscosity coefficient (€) 1s parameterized by a 
dimensionless and adjustable coefficient (N) from equation (21). A first approach for 
the inclusion of turbulent momentum exchange was to determine the minimum least 
square error for spatially-invariant cy and N across a reasonable range of cand N values. 
Bed shear stress coefficients were varied from 0.0005 to 0.0200 and N was varied from 
0.000 to 0.016. A three-dimensional depiction of their solution surface is illustrated in 
Figure 33 for a typical data run. Minimum least square error is represented by a peak 
in the solution surface closest to the top of the three-dimensional box. Note that there 
are many combinations of c, and N that minimize the least square error to nearly the 
same degree. Thus too many degrees of freedom result in multiple mathematical sol- 


utions with a corresponding loss of insight for the physics involved. 


119 


Lf) — @ 
(NJ NN i) en we 
i Ge) Fe CO. eS = 
a ae = = Ve) = 
(oe) Cc) " ) = | 
= O 
II eres ag lI C= ee 
cay eee) « O ee 
C 
uy) 
ee Gx) 
= 
p< 
eo C) 
a Te) 
Lie CO) 
<a 4 
eS ae oS) i ee) 
O I| . | 
ee es es Se 
3) 
uy) 
Lf) ce Wy) a) LY ee Lip C 


Figure 32. Mlodel-generated V compared to field observations versus offshore dis- 
tance -for constant cy, no turbulent momentum exchange, 15-16 
October. Model results are solid line and field observations are “x's. 


Error is mean least square error. 
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The spatially-invariant cy approach is therefore considered a baseline upon 
which to use turbulent momentum exchange to improve the model results. The proce- 
dures are to first determine an invariant c,, and then to vary N until the least square error 
is minimized. This results in a substantial improvement in model results over previous 
approaches. The 18 October data run required an additional c; iteration after a value 
of N was determined. The model result prior to this additional iteration is represented 
by a dashed line in Figure 52. With longshore momentum being transferred over the 
bar and into the trough, the previous velocity underestimation in the trough has been 
significantly reduced. Velocities are modeled within 20% of measured values, with the 
exception of the trough region on 17 October. This case 1s unique and is addressed be- 
low. Interval bars in Figure 31 and Appendix A (Figure 50 through Figure 52) repre- 
sent 10% and 20% of the observed value and are provided solely as a visual aid for 
interpreting the data. Mean least square error ranges between 7-15% forsaliaioes 
transects. 

Current velocities in the trough (Runs 17-6, 17-7, and 17-8) on 17 October are 
much lower (and actually of opposite direction to the longshore current offshore the bar) 
than on the other three days=Guza et al. (1986) showedsinam ims Sy and noteplne 
structure of S,.(/), that is important for longshore current response. It is therefore pos- 
sible that the low current velocities in the trough on 17 October are the result of bi- 
directional wave forcing where breaking swell and sea are driving longshore currents in 
opposing directions and cancelling each other out. Radiation stress spectra (Appendix 
F) show weak wave forcing for the trough region on !7 October. Radiation stress 
spectra for runs 17-1 and 17-2 show sea and swell approaching from different directions 
with relatively large incident angles (18° and —9°, respectivelv). The cross-over wave 
period separating swell from sea is 6.7 seconds. A scatter plot of breaking wave heights 
versus breaking wave periods for a trough run (17-6) shows that approximately half of 
the breaking waves consisted of swell and half consisted of sea (Figure 34). A scatter 
plot for Run 17-4 (located offshore the bar) does not show as much breaking wave ac- 
tivity for the swell. Thus the mintmum combined velocity in the trough on 17 October 
and the poor model results for that particular day and location are possibly a result of 
opposing longshore currents. 

3. Bed Shear Stress Coefficients 
Spatially-invariant c/s are determined from model-fitting using a non-linear 


formulation of ey (TGS86) and turbulent momentum exchange (Table 20). The ratio of 
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Figure 34. Scatter diagram of breaking wave periods and heights for 17 


October. Run 17-4 is offshore the bar and 17-6 is in the trough. 


Table 20. MODEL-FITTING RESULTS FOR B, cp AND N 


ere fon@at 


0.74 0.78 

0.006 OOS 
0.146 O07 
0.003 0.003 
0.003 0.002 


Error is mean least square error for the transect 


‘transect mean cy determined by model-fitting 


2 : 
transect mean cy determined by momentum balance 





wind force to wave force for the SUPERDUCK data was 0.09 (see Chapter VII), there- 
fore the inclusion of surface wind stress in this model did not have significant impact on 
the magnitude of c, from model-fitting for this particular data set. However, as discussed 
earlier, the inclusion of surface wind stress under certain wind, wave, and beach slope 
conditions can have significant impact. Bed shear stress coefficient values for 15-18 
October are 0.003, 0.003, 0.009, and 0.003, respectively. N coefficient values are 0.006, 
0.015, 0.004, and 0.016, respectively. These cy values are in good agreement (sce 
Table 20) with the mean cys determined for each day by the momentum balance 
method. This agreement is encouraging and lends credibility to the cy results considering 
the complex physics being addressed and the different and independent methods em- 


ploved. 
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IX. SUMMARY AND CONCLUSIONS 


Previous investigators of longshore currents have used manv simplifying assump- 
tions which often preclude meaningful and realistic results. The assumptions of a planar 
beach, a steady and depth uniform flow, invariant bed shear stress, and negligible wind 
stress are quantitatively examined by measuring the relative importance of each term in 
a local longshore momentum balance based on measurements acquired during a field 
experiment at a barred beach. Wind and wave forcing of longshore currents is investi- 
gated with a numerical model and from field measurements. Measurements of pressure, 
current velocities at three levels, wind speed, and sea surface slope were acquired froin 
an instrumented sled positioned at various locations across a transect of a barred beach. 
Unfortunately the slope measurements were found during data analysis to be unreliable 
and were disregarded. 

Three-dimensional depictions of the current structure over a nearshore bar are de- 
veloped. The cross-shore velocities [ O(0.2ms~') ] indicate a relative onshore flow in the 
upper third of the water column and a relative offshore flow in the lower third of the 
water column, irrespective of bar location and in agreement with cross-shore theory. 
The longshore velocities [O(1.0ms—') ] indicate nearly depth uniform flow with a slight 
velocity increase with elevation. This depth uniform flow confirms the appropriateness 
of the depth-integration aspect of radiation stress for a barred beach. The longshore 
Vertical current structure did not exhibit a logarithmic profile. 

A numerical model is developed to compare wind and wave forcing for varying wind 
conditions and beach slope. The finite-depth TMA spectral wave model is used as input 
for the energy spectrum. Waves are assumed to be locally generated and fully arisen. 
Wind force is proportional to surf zone width and therefore inversely proportional to 
beach slope. Wave force 1s independent of surf zone width and beach slope. Wind force 
is found to increase more rapidly than wave force for an increase in wind speed or an 
increase in the wind and wave incident angle relative to beach normal. Ratios of wind 
to wave force are calculated. For a given wind direction, a wind speed increase of 10 ms~! 
results in a 50% increase in this ratio. In other words, wind force becomes increasing 
significant relative to wave force, as wind speed and incident wind direction increase. 
For obliquely onshore winds and a steep beach slope (tan Bf = 0.10), wind force is gen- 


erally less than 10% of the wave force and thus may be considered negligible. For ob- 


liquely onshore winds and a gentle beach slope (tan f = 0.01), wind force can be 
sienificant, with wind to wave force ratios ranging from 15% to 100% for winds between 
10-30 ms! and wind and wave directions from 10 — 70°. Thus under certain wind and 
wave conditions, wind force can be a first order term in the longshore momentum bal- 
ance. For offshore winds, longshore currents are fozally wind-driven. In summary, the 
model shows that wind force should be considered as a forcing term for longshore cur- 
rents and that it 1s inappropriate to assume that wind force can be omitted. 

There are no field studies addressing surface wind stress for robust surf zones (H,,, 
greater than 0.5 m). Surf zone atmospheric drag coefficients are determined from surface 
Wind stress measurements acquired immediately outside the surf zone and wind speed 
measurements acquired bv the sled anemometer inside the surf zone. Stability- 
dependent drag coefficients are calculated using the Businger et al. (1971) iterative ap- 
proach. A conservative increase of 33% for C, was made to account for the significantly 
increased surface roughness of the surf zone due to breaking waves based on the results 
of Davidson et al. (1988) which indicate an increased C, for chaotic seas. During periods 
of atmospheric convective mstability and reduced surface roughness (low wave heights), 
Hsu's (1970) empirical relationship for C, is used. Resultant C,s range fron) aaa 
to 3.4x 10-3 for non-convective instability and from 4.6 x 10 to 8:23 xqiG@=aeaee 
convective instability. These C, magnitudes are consistent with previous field measure- 
ments of C, over the open ocean as Well as the theoretical C, values of Geernaert et al. 
(1987) for shallow water depths. 

All terms in the longshore momentuin balance are calculated from field measure- 
ments of pressure, current velocity, and wind speed. The sled was initially positioned 
beyond the surf zone for a complete data run. This additional data run allowed a daily 
numerical calibration of the current meters during data analysis based on conservation 
of éS,,/éx outside the surf zone. This procedure is essential for obtaining realistic 
physical results and is a distinguishing feature of the research. Longshore momentum 
due to unsteady flow was determined to be less than 10% of the longshore momentum 
due to steady flow for this data set. Relative contributions of the radiation stress gra- 
dient (¢S,, / €x), temporal term (¢.1/,/¢1), and surface wind stress (tJ) to the longshore 
momentum balance were calculated and analyzed by location relative to the bar. Mean 
relative contributions for the wave, wind, and temporal terms were 0.86, 0.08, and 0.06, 
respectively, for both offshore and on top the bar. However in the trough, the mean 


relative contributions were significantly different (0.64, 0.17, and 0.19. respectively). For 





the trough region, the relative importance of wave forcing is reduced by 25%, whereas 
the relative importance of surface wind stress 1s doubled and the temporal term relative 
importance is tripled. All three terms are of first order in the trough. Although surface 
wind stress Was nearly constant across the surf zone, the relative contribution of surface 
wind stress increases shoreward from the location of maximum breaking due to the de- 
crease in wave forcing. The relative contribution of the temporal term appears to be 
largest in the trough, possibly due to the presence of surf beat or edge waves. Wind to 
wave force ratios calculated from the field data (0.09) are in agreement with those de- 
termined from the model (0.12) for the same wind conditions. 

Spatially-dependent bed shear stress coefficients are determined between sled run 
locations. The spatial step size (Ax of the radiation stress gradient) 1s determined by 
error analysis to be too small when assigned the 2 m distance between the current meters 
on the sled. It had originally been planned that the slope array, which did not work 
properly, would have provided the necessary precision to measure the radiation stress 
gradient over the 2 m distance. Reasonable c; results are calculated when the spatial step 
size 1s increased to the distance between run positions O(20 - 40 m). Resultant c/s range 
from less than 0.001 to 0.006. Spaually-variable mean c;s are found with values of 0.004 
+ 0.0010 for offshore the bar, 0.002 + 0.0006 for on top the bar, and 0.001 + 0.0003 for 
the trough (excluding the nearshore slope). A trend for a shoreward decreasing Cc, 1s 
noted, possibly due to the shoreward decrease in wave action. A mean cy for the 
nearshore slope is not specified since error analvsis shows that a significant overesti- 
mation of c, can occur for the nearshore slope due to the combined effect of decreased 
depth and steeper slope upon the fimite difference approximation. The percentage 
change in cy due to neglecting surface wind stress of the temporal term in the momentum 
balance ranges from near zero to 100%, with a mean change between 10-30%. Thus for 
this data set and under certain conditions, the surface wind stress term and the temporal 
term can be of first order in the momentum balance. 

Longshore current velocities at a barred beach are modeled to within + 20% of the 
observed values at SUPERDUCK bv modifying the Thornton and Guza (1986) model 
to include wind stress, wave incident angle dependency, slope-dependent y, and an im- 
proved breaking wave weighting function. Mean least square error between modeled 
and observed current velocity ranges from 7-15%. The inclusion of turbulent momen- 
tum exchange to transport momentum into the more placid trough region is required 


and is physically reasonable. Spatially-invariant cy values determined by model-fitting 


ey 


for 15-18 October are 0.003, 0.003, 0.009, and 0.003, respectively. These cy Values agree 
well with mean cy values determined for each day by the momentum balance method. 
This agreement is encouraging and lends credibility to the cy results considering the 


complex physics being addressed and the different and independent methods employed. 
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APPENDIX A. DATA ANALYSIS FIGURES FOR DAYS OTHER THAN 
“OCTOBER 16, 1986 
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Figure 35. Mean current structure on 11 October 1986. Numbers represent data 


runs, short dotted lines are MWL’s, and distance alongshore is arbi- 
trary. 
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Figure 36. Nfean current structure on 12 October 1986. Numbers represent data 


runs, short dotted lines are MWL’s, and distance alongshore is arbi- 
trary. 
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Figure 37. MNlean current structure on 15 October 1986. Numbers represent data 
runs, short dotted lines are MWL’s, and distance alongshore is arbi- 


trary. 
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Figure 38. Nfean current structure on 17 October 1986. Numbers represent data 


runs, short dotted lines are MWL’s, and distance alongshore is arbi- 
trary. 
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Figure 39, Nilean current structure on 18 October 1986. Numbers represent data 


runs, short dotted lines are MWL’s, and distance alongshore 1s arbi- 
trary. 
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Figure 40. Mfean current structure on 21 October 1986. Numbers represent data 


runs, short dotted lines are MWL’s, and distance alongshore 1s arbi- 
trary. 
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Figure 41. Offshore S,, spectra and S,J from ‘puw’ gage 621 for 15 October 


1986. Numbers indicate data runs. Intervals are 95% confidence. 
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Figure 42. Offshore S,, spectra and SJ from ‘puy’ gage 621 for 17 October 
1986. Numbers indicate data runs. Intervals are 95% confidence. 
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Figure 43. Offshore S,, spectra and S,/ from ‘puv’ gage 621 for 18 October 
1986. Numbers indicate data runs. Intervals are 95% confidence. 
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Figure 44. 
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Wind speed, wave height, atmospheric stability, and initial C, data for 
October 11-12, 1986. Beach is assumed parallel to the y-axis. Open 
arrowheads denote wind direction and closed arrowheads denote wave 


incident angle. Winds are adjusted to z = 10 m. 
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Wind speed, wave height, atmospheric stability, and initial C, data for 
October 17-18, 1986. Beach is assumed parallel to the y-axis. Open 
arrowheads denote wind direction and closed arrowheads denote wave 


incident angle. Winds are adjusted to z = 10 m. 
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Figure 46. Wind speed. wave height. atmospheric stability, and initial C, data for 
October 21, 1986. Beach is assumed parallel to the y-axis. Open 
arrowheads denote wind direction and closed arrowheads denote wave 


incident angle. Winds are adjusted to z = 10m. 
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Figure 47. Wave height and breaking wave height distributions (hatched area) nor- 


malized to offshore H/,,. for 15 October... The Rayleigh distribution 1s 


rms 
the solid line and empirical breaking wave distributions given by 


equations (127) and (140) are the dashed and dotted lines, respectively. 
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Figure 48. Wave height and breaking wave height distributions (hatched area) nor- 
malized to offshore H,,,, for 17 October. The Ravleigh distribution 1s 
the solid line and empirical breaking wave distubutions given by 
equations (127) and (140) are the dashed and dotted lines, respectively. 
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Figure 49. Wave height and breaking wave height distributions (hatched area) nor- 
malized to offshore H,,,, for 18 October. The Ravleigh distribution is 
the solid line and empirical breaking wave distributions given by 
equations (127) and (140) are the dashed and dotted lines, respectively. 
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Figure 50. Model-generated H,,,, and V compared with field observations versus 


distance for invariant cy and B and including turbulent momentum ex- 
change, 15 October. Model results are solid lines and field observa- 


tions are “x”’s. Elevation is referenced to MSL. 
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Figure 51. 





Nlodel-generated 


rms 


and V compared with field observations versus 


distance for invariant cy and B and including turbulent momentum ex- 
change, 17 October. Model results are solid lines and field observa- 


tions are “x’’s. Elevation is referenced to MSL. 
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Model-generated H,,,. and V compared with field observations versus 
distance for invariant cp and B and including turbulent momentum ex- 
change, 18 October. Model results are solid lines and field: observa- 
tions are “x”’s. Elevation is referenced to MSL. Model result prior to 


the last c, iteration is represented by the dashed line. 
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Figure 53. Nfodel-generated V compared to field observations versus offshore dis- 
tance for constant cy and no turbulent momentum exchange, 17-18 Oc- 
tober. Model results are solid line and field observations are “x”’s. 


Error 1s mean least square error. 
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APPENDIX B. CURRENT VECTORS OVERLAYED ON BATHYMETRY 
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Sled and alongshore array mean current vectors for 11 October 1986 


Figure 54. 


Velocity scale is at upper right. 


overlayed on bathymetry. 
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Sled and alongshore array mean current vectors for 12 October 1986 


Figure 55. 
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Velocity scale is at upper r1 


overlaved on bathvmetry. 
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Figure 


Sled and alongeshore array mean current vectors for | 


Velocity scale is at upper right. 


overlaved on bathymetry. 
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Sled and alongshore array mean current vectors for 16 October 1986 


Figure 57. 


Velocity scale is at upper night. 


overlaved on bathymetry. 


153 


eS 
— 


Use : 
(HH) GYOHSSNOW JONYULSIO 


USUI 


US6 





DISTANCE Ghee Segre tine 


Sled and alongshore array mean current vectors for 17 October 1986 
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Figure 


Velocity scale is at upper right. 


overlaved on bathymetry. 
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Sled and alongshore array mean current vectors for 18 October 1986 
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Figure 


Velocity scale is at upper right. 


overlayed on bathymetry. 
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Sled and alongshore array mean current vectors for 21 October 1986 


Figure 60. 


Velocity scale is at upper right. 


overlaved on bathymetry. 
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APPENDIX C. ATMOSPHERIC STABILITY DETERMINATION 


NiecwiolliGmang Set Of Equations represent an iterative procedure to determine z,L 


based on the stability flux profile relation of Businger et al. (1971). 


Initial equations: 


Tew = Tair + 35( In(.01RH)) (C1) 
AT = Top — Tyurface — -0098Z7 (C2) 
GG = Ree eR bs (C3) 


Assign all Tae = it) 


Deternune 7” from iteration: 


ne (Ca) 





Ya= Fo | Y= 2 Inf Base) 5 0 (C2) 
Y =-7(+) =o 
L J 
O =O Trew + 273.15)(.01 RA) (C6) 
OO 27215) (C7) 
unalone 23 84 5 03 log(T..) (C8) 
Kel 

* (Q ate Qo)K 

eS (C9) 
In Zo Ea a 
i eo leah0n (iO ) (C10) 


Sy 


= (C11) 


SF Fe seal Ds (C12) 
P &K Zo ie 

7 a ar (Gig) 
Kel 


Continue iterations from equation (C4) until Ca) converges. 


Definitions: 


temperature measurement height (m), 
humidity measurement height (1m), 


relay MUmidiove 6), 


= dewpoint temperature (degrees C) 


= temperature roughness length (2 x 107°) 


humidity roughness length (2 x 107°) 


= integral diabatic term for temperature 


integral diabatic term for humidity 


specific humidity (g kg) 


APPENDIX D. SURFACE WEATHER MAPS 
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Figure 61]. 
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for 0700 E.S.T. 16 October 1986. 


Surface weather analysis 


Figure 62. 
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for 0700 E.S.T. 17 October 1986. 


Surface weather analysis 


Figure 63. 


162 








701 


Bie 


821. 





~— - 


pos! 


ere, Aeon PB * ~~) 


Comin 


“Fb he DVM Brtaeee ga Pa bedy e 


OF PON. 10% 00 BAe, 80.0270 toe 


Z101 


bh 
eee Ny : 
-. 


ae 1 ve 








18 October 1986. 


He 


S 


ysis for 0700 E. 


Surface weather analysis 


Figure 64. 
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APPENDIX E. ANALYSIS PARAMETERS FOR ALL SLED DATA RUNS 
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® x,y coordinates are referenced to the FRF coordinate system 





° f/f, \8 the peak freq m the sea sfc elevation spectrum 


© “EMW’' is the “equivalent monochromatic wave” incident angle (Thornton and 
Guza, 1986) 


e S,,/ is the mean S,,/ for the submerged current meters on that run 


Jm 1S the peak frequency of the sea surface elevation spectra 
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APPENDIX F. RADIATION STRESS, SEA SURFACE ELEVATION, 
SIGNIFICANT WAVE ANGLE, AND CURRENT VELOCITY SPECTRA 
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Figure 65. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(z), and current velocity spectra for sled run 15-1. 
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Figure 66. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


A ° - 
(x), and current velocity spectra for sled run 15-2. 
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Figure 67. Radiation stress (S,.), sea surface elevation (7), significant wave angle 


(x), and current velocity spectra for sled run 15-3. 
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Figure 68. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(z), and current velocity spectra for sled run 15-4. 
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Figure 69. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(x). and current velocity Spectra for sled run 16-2. 
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Figure 70. 


Radiation stress (S,,), sea surface elevation (7), significant wave angle 
(2), and current velocity spectra for sled run 16-3. 
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Figure 71. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(x), and current velocity spectra for sled run 16-4. 
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Figure 72. Radiation stress (S,.), sea surface elevation (7), significant wave angle 


(x), and current velocity spectra for sled run 16-5. 


al 


RUN 16-6 —-—- SUPERDUCK 
O57 y (OV ons) 


yx a 
3500 
~-en 
re 
'e)) 
= 2 8 : 
/) 
<= ai 
> a 
<< 500-4 Su | 
© : ene ere een GR 1 Oo 
DEOO.1° 0.2055 0240's Hi. 1 0.2 0.3 O145Gre 
1000-1 al 


G(f) (m*s~") 





0 = eee 
0.00.1 0.2 0.3 0.4 0.5 
Frequency (Hz) 


| | \ 
0.00.1 0.2 0.3 0.4 0.5 
Frequency (Hz) 


Figure 73. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(z), and current velocity spectra for sled run 16-6. 
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Radiation stress (S,x), sea surface elevation (17), significant wave angle 


(z), and current velocity spectra for sled run 17-1. 
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Figure 75. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(z), and current velocity spectra for sled run 17-2. 
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Figure 76, Radiation stress (S,.,), sea surface elevation (7), significant wave angle 


(zx), and current velocity spectra for sled run 17-3. 
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Figure 77. 


Radiation stress (S,,), sea surface elevation (7), significant wave angle 
(x), and current velocity spectra for sled run 17-4. 
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Figure 78. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(a), and current velocity spectra for sled run 17-5. 


183 


RUN 17-6 ——" SUPERP ee. 
05% Xx (One) 








Dyx OX 
| £55) 6(6) 30 
ca 
oo, 
| 2000 S28 
< <b) 
7% a 
> 
Aaa 
< 500 = 
oe @S 
OC (= an 
| 0°50 0: Le Oerne sO) Coes 
-100-! ~10 
T) Current Velocity 
1 Vel 
U 
V 





| \ j ’ 
0.00.1 O2296:5 10.4485 


! ae ‘ac 
0.00.1 0.2 0.3 0.4 0.5 
Frequency (Hz) 


Frequency (Hz) 


Figure 79, Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(x), and current velocity spectra for sled run 17-6. 
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Figure 80. 


Radiation stress (S,,), sea surface elevation (7), significant wave angle 
(z), and current velocity spectra for sled run 17-7. 


185 


G(f) sm“) 


G(f) (m’s) 


Figure 81. 


RUN 17-8 —— SUPaRewien 
95% x (O74 








Dyx OX 
S00 35 
Soa 
2000 D) 20 
<2 
] 10 
LS 
> Sa 
500- 5 : A AN 
i { i : ; } 0 k | wv i Fin fi 
0.00.1 0.2 0.3 0.4 0.5 00.1 0,2 O15 Oris 
-1000- kaa 
T) Current Velocity 
3 Vel 
U 
V 


G(f) (m’s7') 





oe ar ‘ a: vk an i iT 
0.00.1 0.2 0.3 0.4 0.5 0.00.1 0.2 0.3 0.4 0.5 
Frequency (Hz) Frequency (Hz) 


Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(a), and current velocity spectra for sled run 17-8. 
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Figure 82. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(z), and current velocity spectra for sled run 18-1. 
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Figure 83. 


Radiation stress (S,,), sea surface elevation (7), significant wave angle 
(x), and current velocity spectra for sled run 18-2. 
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Figure 84. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(z), and current velocity spectra for sled run 18-3. 
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Radiation stress (S,,), sea surface elevation (7), significant wave angle 
(z), and current velocity spectra for sled run 18-4. 
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Figure 86. Radiation stress (S,,), sea surface elevation (7), significant wave angle 


(x), and current velocity spectra for sled run 18-5. 
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APPENDIX G. MEAN CURRENT VELOCITIES AND TOTAL 
RADIATION STRESS FOR EACH SLED CURRENT METER AND RUN 
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